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SYSTEM  ThOg-YgO^* 


Eo  Co  Subbarao  and  F.  Ho  Sutter 


ABSTRACT 


Compositions  in  the  ThOg-YO^  ^  system  were  eo-precipitated 
as  the  oxalates.  Disks,  pressed  from  decomposed  oxalates,  were 
sintered  in  oxygen  at  1400-2200  C.  Densities  of  the  sintered  disks 
were  96-98#  of  the  theoretical.  Solid  solutions  with  the  fluorite- 
type  structure  are  formed  up  to  20-25  mole  #  YO.,  j-  at  1400  C  and  up 
to  45-50  mole  #  YGLj  ^  at  2200  C.  Density  da tajj have-shown-  that  these 

'  t  1 

solid  solutions  correspond  to  Tbi_xYx02_o  havin^  8  complete 
cation  sublattice .  filled  by  the*  andA  Y^ions,  and  vacancies 

in  the  anion  sublattice .  The  observed  increase  in  electrical 
conductivity  with  increase  in  YO^  ^  content  is  consistent  with  the 
notion  that  oxygen  ions  are  the  predominant  charge  carriers  through 
a  vacancy  mechanism.  . 


JvS— 0ff±ce-of-Nav3l-Research._ 
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INTRODUCTION 


Thorium  dioxide,  ThOg,  has  the  cubic,  CaFg-type  structure  up 

to  its  melting  point,  ~330O°C.  yttrium  sesqui oxide,  YgO^  (or  Y01  ^), 

crystallizes  with  the  C-type  rare  earth  oxide  structure  (or  TlgO^  type). 

This  structure  is  closely  related  to  that  of  CaFg,  from  which  it  may 

be  derived  by  an  ordered  omission  of  one-quarter  of  the  anions  and 

then  a  slight  rearrangement  of  the  atoms.  Based  on  the  similarity  of 

the  structure,  extensive  solid  solution  may  be  expected  in  the  binary 

system  ThO„-YO,  •  Hund  and  Mezger^  found;  that  22  to  30  mole  $  YO,  - 
d.  JL*p  J-O 

is  soluble  in  ThOg  at  1200°C. 

The  purpose  of  the  present  investigation  is  to  determine  the 
phase  relations  and  solid  solubility  in  the  binary  system  Th0g-Y01  ^  at 
temperatures  up  to  2200°C,  to  examine  the  defect  structure  of  the 
resulting  solid  solutions  and  to  present  a  preliminary  account  of  the 
correlation  between  the  defect  structure  and  ionic  conductivity  measured 
on  dense  ceramic  disks.  Further  details  of  the  ionic  conductivity  of 
ThOg  solid  solutions  will  be  discussed  elsewhere. 

II.  SPECIMEN  PREPARATION 

Thorium  nitrate  and  yttrium  nitrate,  obtained  from  Lindsay 
Division  of  the  American  Potash  and  Chemical  Corporation,  are  the 
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starting  materials.  Spectrographic  analyses  of  these  materials  are  given 
in  Table  I.  The  nitrates  were  dissolved  in  distilled  water  and  the  oxide 
contents  were  determined  by  ignition  at  1300-l400°C.  The  nitrate  solutions 
were  mixed  in  the  required  proportions  to  prepare  the  desired  solid  solu¬ 
tion  compositions.  Thorium  and  yttrium  salts  can  be  precipitated  either 
as  (A)  the  hydroxide  or  (B)  the  oxalate. 

(A)  Hydroxide:  In  this  method,  described  in  detail  by  Hund  and 
Mezger, ^  the  heated  nitrate  solution  was  added  to  boiling  ammonium 
hydroxide  solution  with  vigorous  stirring.  The  precipitate  was  filtered, 
washed  and  dehydrated  at  105°C.  After  grinding,  it  was  decomposed  to 

the  oxides  at  1200°C.  The  oxide  mixture  was  ground,  pressed  into  disks 
and  the  disks  were  sintered  at  temperatures  between  l4O0°  and  2000°C.  The 
densities  of  disks  that  were  prepared  from  the  hydroxide  and  from  the 
.oxalate  are  compared  in  Table  II.  The  densities  were  determined  from 
the  mass  and  the  dimensions.  Inasmuch  as  the  oxalate  method  led  to 
higher  densities,  the  hydroxide  method  was  discarded. 

(B)  Oxalate:  Since  the  particle  size  distribution  of  the 

starting  oxide  powder  has  a  remarkable  influence  on  the  sintering 

behavior  of  the  oxide,  it  is  advantageous  to  obtain  fine-grained  oxides. 

In  a  sintering  study  of  thoria  obtained  from  the  hydroxide,  carbonate, 

(2) 

nitrate,  and  oxalate,  Kantan  et  al  concluded  that  the  oxalate  method 
yielded  the  maximum  density. 


I 


-  k  - 

In  this  method.,  the  nitrate  solutions  were  placed  in  an  ice 
hath.  Saturated  oxalic  acid  was  added  drop  by  drop  while  the  chilled 
nitrate  solution  was  vigorously  stirred.  Beyond  this  step,  two  courses 
were  available;  (l)  dehydration  of  the  entire  solution,  or  (2)  filtra¬ 
tion  and  dehydration  of  the  precipitate. 

(1)  The  entire  solution  was  dehydrated  at  about  120°C, 

leaving  a  cake  at  the  bottom  of  the  beaker.  The  top  surface  of  the 

cake  was  more  granular  than  the  bottom.  This  difference  in  grain  size 

became  more  pronounced  on  heating  at  about  900°C.  X-ray  diffraction 

studies  of  the  different  portions  of  the  cake  indicated  the  granular 

material  at  the  top  of  the  cake  to  be  mostly  YO,  c  and  the  fine  material 

-*•0 

at  the  bottom  to  be  essentially  ThOg.  The  reasons  for  segregation  and 

grain  growth  of  Y01  are  not  clear  at  present.  Disks,  pressed  from 

x.p 

this  powder  without  much  grinding,  had  an  inhomogeneous  appearance 
after  firing  to  2000°C  for  4  hrs.  White  specks  (of  YgO^)  were  randomly 
distributed  throughout  the  volume  of  the  disk.  Careful  grinding  was 
necessary  to  obtain  homogeneous  solid  solutions  from  material  prepared 
by  this  method. 

(2)  The  precipitate  was  filtered,  taking  care  to  avoid 
segregation  in  the  beaker,  washed  until  the  precipitate  was  free  of 
acids,  and  then  dehydrated  at  about  120°C.  The  resulting  fine-grained 
oxalates  were  calcined  to  the  oxides  at  900°C^^  for  8-12  hrs.  The 


oxide  powder  was  homogeneous  and  fine-grained.  This  method  gave  satis¬ 
factory  results  and  was  therefore  followed  in  the  rest  of  the  investigation. 

The  powder  was  ground  and  pressed  into  disks,  3-12.5  mm  in  diam. 
and  1-10  mm  in  thickness  at  a  pressure  of  about  20,000  psi.  After  pre¬ 
calcination  at  1400°C  for  4  hrs.  in  a  globar  furnace,  the  specimens  were 
sintered  in  an  oxygen  atmosphere  in  the  temperature  range  1600  to  2200°C 
for  3  to  7  hrs.  using  a  propane-oxygen  pot  furnace.  During  sintering, 
the  specimens  were  enclosed  in  covered  thoria  crucibles  to  avoid  contact 
with  undesirable  vapor  species  that  may  be  present  in  the  furnace  chamber. 

To  prevent  the  disks  from  sticking  to  each  other  during  sintering,  they 
were  separated  using  powders  sintered  at  2250°C  as  spacer  material. 

The  ThOg  solid  solution  compositions  sintered  in  oxygen  range 
in  color  from  ivory  white  to  brown,  the  intensity  of  the  brown  color  in¬ 
creasing  with  YiQ^  j.  for  a  given  sintering  temperature  and  increasing 
with  temperature  for  a  given  TO ,  content. 

1.5  - 

Typical  microstructures  of  polished  and  etched  surfaces  of 
the  sintered  specimens  are  shown  in  Fig.  1.  The  pure  thoria  was  etched 
by  hot  phosphoric  acid  and  the  others  by  hot  sulfuric  acid.  The  grain 
size  of  all  of  the  specimens  is  of  the  order  of  20  microns,  and  increases 
slightly  with  increasing  YD\  _  content.  The  pore  size  also  appears  to 
increase  with  YO^ . ^  content.  The  porosity  has  been  roughly  estimated 
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at  2-1*  percent  in  the  specimens  with  0  to  15  mole  percent  YC1  the 

higher  porosity  applying  to  the  higher  YO,  c  content. 

■*■0 

III.  PHASE  RELATIONS 

Specimens  that  had  been  sintered  between  ll*00°  and  2200°C 
for  3  to  6  hrs.  in  an  oxygen  atmosphere  were  cooled  in  the  furnace 
after  turning  off  the  gas  flame.  X-xay  examination  of  the  sintered 
specimens,  using  Cu  Ka  radiation,  was  used  to  determine  the  phases 
present. 

The  system  may  conveniently  be  divided  into  three  regions: 
ThOg-type  solid  solutions,  YgO^-type  solid  solutions  and  a  two  phase 
region  separating  these  two  solid  solutions.  Solid  solutions  with 
fluorite-type  structure  were  found  up  to  1*5-50  mole  £  Y0  ^  -  at  2200°C, 
decreasing  to  20-25  mole  %  YO,  c  at  ll*00°C.  The  latter  value  is 
consistent  with  Hund  and  Mezger’s  results. ^  The  lattice  parameters 
of  specimens  sintered  at  2200°C  are  given  in  Table  III  and  plotted  in 
Fig.  2.  These  values  were  obtained  from  diffraction  lines  with  20  >  100° 
and  extrapolated  according  to  the  method  of  Nelson  and  Riley. ^  The 
decrease  of  lattice  constant  with  YO,  c  content  is  due  to  the  smaller 
size  of  the  Y^*  ion  (0.92  A)  compared  to  the  Th^  ion  (0.99A). 
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IV.  DETECT  STRUCTURE  AND  IONIC  CONDUCTIVITY 
OP  THE  FLUORITE  PHASE 

4+ 

When  the  Th  ions  are  replaced  by  Y  ions,  charge  neutrality 

can  be  restored  by  one  of  two  structural  schemes:  (l)  The  cation  sub- 

4+  3+ 

lattice  is  completely  filled  by  Th  and  Y  ions  and  the  appropriate 
number  of  anion  vacancies  are  created.  (2)  The  anion  sublattice  is 
complete  and  the  excess  cations  occupy  interstitial  sites.  In  order 
to  distinguish  between  these  two  models,  densities  corresponding  to 
them  are  computed  using  the  lattice  parameters  listed  in  Table  III. 

These  density  values  are  compared  in  Fig.  3  with  those  measured 
pycnometrically  on  powders,  sintered  at  2200°C,  using  butyl  phthalate 
as  the  liquid.  On  the  basis  of  these  data,  it  was  concluded  the  ThOg 
solid  solutions  have  fluorite  type  structure  with  a  filled  cation 
sublattice  and  anion  vacancies.  These  solid  solutions  may  therefore 

be  represented  by  Th^Va-0. 5x’ 

The  densities  of  sintered  disks  used  for  conductivity 
measurements  were  measured  by  a  displacement  method  using  butyl 
phthalate.  These  values  are  also  included  in  Fig.  3  and  show  that 
these  specimens  have  a  porosity  of  1-%,  in  agreement  with  the  photo¬ 
micrographic  study. 

Electrical  conductivity  was  measured  on  the  sintered 
disks  with  platinum  paste  electrodes  at  a  frequency  of  1000  cps 


vising  a  modified  General  Radio  capacitance  bridge  in  conjunction  with 
decade  resistance  and  capacitance  boxes.  The  data  for  ThO-  (obtained 
by  the  decomposition  of  the  nitrate  as  well  as  the  oxalate)  and  for 
ThOg  with  1,  2  and  %  YO^  for  the  temperature  range  500°  to  l4oo°C 
are  shown  in  Fig.  U,  where  log  oT  is  plotted  versus  reciprocal  absolute 
temperature.  The  conductivity  increases  with  increasing  YO,  content. 

-L.  p 

This  behavior  is  consistent  with  the  notion  that  the  current  carrying 
species  in  these  solid  solutions  are  the  anion  vacancies  and  that  the 
conductivity  increases  with  anion  vacancy  concentration. 
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Table  I  —  Spectrographic  Analyses  of  Starting  Materials 


TH(No3)4  *^03)3 


Element 

PPP 

Element 

PPm 

Bi 

<  1 

Mg 

20 

Ca 

2.5 

Si 

200 

Cu 

2 

Pe 

200 

In 

<  1 

Mn 

50 

Mg 

3 

Ca 

<  100 

Nb 

<  2 

A1 

50 

Fb 

<  1 

Cu 

<  100 

Si 

15 

Ti  - 

50 

Sn 

<  1 

Co 

<  100 

Ti 

<  2 

Ni 

<  100 

U 

<  100 

W 

<  30 

Zn 

<  10 

Zn 

3 

Y 

<  2 

Ce 

<  200 

Au 

<  10 

Ba 

<  3 

K 

<  6 

Na 

3 

P 

<  50 

Sb 

<  1 

Sn 

<  1 
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Table  II  —  Densities  of  Sintered  Thoria  Specimens 


Sintering 

Conditions 

Hydroxide 
(Calcined  at 

600°C  -  8  hrs. ) 

Oxalate 
(Calcined  at 
900°C  -  12  hrs. ) 

1200°C  -  2  hrs. 

J.k  g  cm"^ 

7.3  g  CD 

lU00°C  -  2  hrs. 

7.5 

7-8 

l600°C  -  2  hrs. 

7-3 

8.9 

1800°C  -  2  hrs. 

7.9 

9-7 
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Table  III  —  Lattice  Parameters 


Mole  $  YQlt5» 


Lattice  Constant,  a,  8 


0 

5-596 

±  0.001 

5 

5-568 

±  0.001 

10 

5.582 

±  0.001 

15 

5.572 

±  0.001 

20 

5.568 

±  0.001 

25.7 

5.549 

±  0.001 

30 

5.540 

±  0.001 

35 

5.527 

±  0.001 

4o 

5.512 

±  0.001 

50 

5.490 

±  0.001 

55 

5-483 

±  0.002 

6o. 

5.500 

±  0.002 

^-Specimens  were  sintered  in  oxygen  at  2200°C  for  3  hrs. 
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FLUORITE  PHASE  IN  THE  SYSTEM  ZrOg-CaO: 
X-RAY  AND  ELECTRICAL  CONDUCTIVITY  STUDY* 


by 


T.  Y.  Tien  and  E.  C.  Subbarao 


ABSTRACT 


In  the  ZrOg-CaO  system  cubic  solid  solutions  of  the  fluorite- type 
structure  are  formed  between  12  and  22  mole  per  cent  CaO  for  specimens 
fired  at  2000°C  and  subsequently  quenched  from  l400°C.  The  lattice  para¬ 
meter  increased  linearly  from  5*125  A  to  5.144  A  with  increasing  CaO 
content.  Results  of  X-ray  diffraction  intensity  studies  favored  the 
oxygen  ion  vacancy  model  over  the  cation  interstitial  model  for  these 

solid  solutions.  As  the  CaO  content  was  increased  from  13  to  20  mole  per 

-2 

cent,  the  electrical  conductivity  at  1000 "c  decreased  from  5  x  10  to 
6  x  10  J  ohm  cm  and  the  activation  energy  increased  from  1.11  to  1.35 
e.v.  Specimens  heated  at  1000 °C  for  one  week  showed  extra  X-ray  diffraction 
lines.  These  data  coupled  with  accompanying  changes 'cin  electrical  conduc¬ 
tivity  are  discussed  in  terms  of  an  order- disorder  transition. 


*  Supported  by  ARPA  through  Office  of  Naval  Research 
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I ,  INTRODUCTION 

Over  the  years,  zirconia,  ZrO^,  and  zirconia  combined  with  oxides 
of  di-  and  tri-valent  metals  have  been  the  subject  of  numerous  investi¬ 
gations,  first  as  a  refractory  and  more  recently  as  a  solid  electrolyte 
1  2 

in  galvanic  and  fuel  cells.  Summaries  of  the  voluminous  literature  on 

3  4  5 

zirconia  are  presented  by  Ryshkewitch,  Roth,  Stocker,  and  Collongues 
et  al.^  among  others. 

Zirconia  has  monoclinic  symmetry  at  room  temperature  and  transforms 

to  a  tetragonal  form  at  about  1100 °C,  accompanied  by  a  large,  disruptive 

7 

volume  change.  Duwez  et  al.  showed  that  zirconia  compositions  contain¬ 
ing  16  to  30  mole  per  cent  CaO  have  cubic  symmetry  when  quenched  from 

Q 

2000 °C.  On  the  otherhand,  Hund  found  the  cubic  phase  to  exist  from  10 
to  20  mole  per  cent  CaO  in  specimens  prepared  at  l460°C.  Dietzel  and 
Tober^  reported  that  cubic  solid  solutions:  pxtend  from  7  to  24  mole  per 
cent  CaO  at  1800°C  and  from  14  to  20  molp  per  cent  CaO  at  l400°C.  A 
recent  Russian  investigation10  places  the  cubic  phase  field  between  10 
and  40  mole  per  cent  CaO  in  specimens  prepared  at  1500°C.  Thus,  consid¬ 
erable  disagreement  exists  about  the  cubic  phase  region  in  the  ZrOg-CaO 
system. 

The  cubic  solid  solutions  in  this  system  have  the  fluorite-type 
structure .  Hund^  compared  measured  densipi.es  with  values  computed  for 
two  structure  models :  one  in  which  anion  vacancies  are  created  and  the 
other  in  which  excess  cations  occupy  interstitial  sites.  His  results 
support  the  anion  vacancy  model,  i.e.  the  Zr  apt  ions  completely 

—fill  the  cation  sites  and  electrical  neutrality  is  preserved  by  the 
creation  of  anion  vacancies,  the  number  of  vacancies  created  being  equal 
to  the  number  of  molecules  of  CaO  added. 

3  11 

Electrical  conductivity  measurements  by  Hund,  Trombe  and  Foex, 

1  O 

Kingery  et  al.  ^  and  others  revealed  rather  high  electrical  conductivities 
for  these  cubic  solid  solutions  at  elevated  temperatures.  Kingery  et  al. 
have  demonstrated  by  means  of  oxygen  diffusion  and  conductivity  measurements 
that  the  transference  number  for  oxygen  ions  in  ZrQ<g,-  Ca^,.  O^g^  is 
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near  unity.  Recently,  Rhodes  and  Carter  measured  the  cation  mobility 

in  these  solid  solutions  and  found  that  the  diffusion  rates  of  both 
4+  2+ 

Zr  and  Ca  ions  were  several  orders  of  magnitude  lower  than  that 

l4 

measured  for  the  oxygen  ion.  Weissbart  and  Ruka's  results  indicate 
that  the  electronic  conductivity  makes  rather  a  small  (perhaps  less  than 
2  per  cent)  contribution  to  the  total  conductivity.  Thus,  the  electri¬ 
cal  conductivity  is  essentially  due  to  the  movement  of  oxygen  ions  through 
anion  vacancies.  This  relationship  between  conductivity  and  anion  vacancies 
in  these,  solid  solutions  might  be  expected  to  lead  to  an  increase  in  conduc¬ 
tivity  with  an  increase  in  the  concentration  of  anion  vacancies  (and  CaO 
content).  However,  the  limited  available  data  such  as  that  reported  by 
Trombe  and  Foex"^  indicate  decreasing  conductivity  with  increasing  CaO 
content  in  the  cubic  phase  field. 

The  present  investigation  was  undertaken:  (l)  to  delineate  the  cubic 
phase  boundaries  when  high  purity  materials  are  used,  (2)  to  obtain  support 
for  the  anion  vacancy  structure  model  for  the  fluorite-type  phase  by  means 
of  X-ray  diffraction  intensity  studies,  (3)  to  measure  electrical  conduc¬ 
tivity  of  several  cubic  solid  solution  compositions  over  a  wide  range  of 
temperature,  and  (4)  to  attempt  an  explanation  for  the  relationship  be¬ 
tween  conductivity  and  anion  vacancy  concentration. 

II.  SPECIMEN  PREPARATION 


Chemically  pure  calcium  carbonate  and  hafnium-free  ZrO^  were  used 
to  prepare  the  CaO-ZrO^  compositions.  The  spectrographic  analysis  of 
the  ZrO^  appears  in  Table  I.  The  weighed  mixtures  were  wet  blended 
using  acetone  as  the  vehicle,  dried,  and  then  reacted  at  1350 °C  for  24 
hours  in  platinum  containers.  The  reacted  mixtures  were  comminuted,  com¬ 
pacted  into  right  circular  cylindrical  specimens  and  then  sintered  in  an 
oxygen  atmosphere  in  a  muffled  propane- oxygen  pot  furnace  at  2000 °C  for 
2  hours.  The  sintered  specimens  were  subsequently  annealed  at  l400°C  for 
one  week.  The  apparent  densities  of  all  compositions  were  in  the  range 
90  to  95  per  cent  of  theoretical.  The  theoretical  densities  were  computed 
from  lattice  parameter  data  and  were  based  on  the  oxygen  vacancy  model. 

The  representative  microstructures  are  presented  in  Figure  1. 


III.  RESULTS  AND  DISCUSSION 


1.  Solid  Solubility 

X-ray  diffraction  photographs  vere  obtained  with  a  11.4  cm  Debye 
-Scherrer  camera  on  a  Norelco  diffraction  unit  using  Cu  kk  radiation. 
Single  phase  cubic  fluorite-type  solid  solutions  were  observed  for  com¬ 
positions  containing  approximately  12  to  22  mole  per  cent  CaO.  These 
specimens  were  heated  at  2000 *C  for  2  hours,  annealed  at  l400°C  for  one 
week  and  then  air  quenched  from  l4O0°C.  The  lattice  parameters  versus 
composition  are  plotted  in  Figure  2.  The  values  were  obtained  by  measur¬ 
ing  diffraction  lines  with  2  6 > 100°  and  extrapolating  by  the  method  of 
Nelson  and  Riley.1'*  The  back  reflection  lines  of  compositions  contain¬ 
ing  less  than  12  mole  per  cent  CaO  were  always  so  diffuse  that  accurate, 
lattice  parameter  measurements  could  not  be  made.  Compositions  contain¬ 
ing  11  mole  per  cent  CaO,  or  less  showed  monoclinic  ZrO^*  CaZrO^  was 
detected  in  samples  with  24  mole  per  cent  CaO  or  more.  Thus  the  X-ray 
study  places  the  limits  of  the  cubic  phase  at  approximately  12  and  22 
mole  per  cent  CaO. 

Crushed  fragments  of  the  specimens  were  examined  with  a 
petrographic  microscope  using  crossed  nicols.  Birefringent  phases  (mono¬ 
clinic  ZrOg  and  CaZrO^)  were  detected  in  12  and  22  mole  per  cent  CaO 
compositions.  Polished  and  etched  surfaces  of  the  sintered  specimens 
were  viewed  under  a  metallographic  microscope.  Traces  of  second  phase 
inclusions  were  observed  in  12  and  22  mole  per  cent  CaO  samples,  as  shown 
in  Figures  1  a  and  1  c. 

According  to  the  present  results,  the  cubic  phase  extends 
from  between  12  and  13  mole  per  cent  to  between  20  and  21  mole  per  cent 
CaO  under  the  conditions  of  preparation  employed.  These  limits  are  in 
fair  agreement  with  Hund's  data. 

2.  X-Ray  Intensity  Study 

The  Zr0o-Ca0  solid  solutions  crystallize  with  the  cubic 
c  8 

fluorite-type  structure  of  space  group  Fm3m.  Hund'  proved  by  density 

measurements  that  the  cation  lattice  sites  were  completely  filled  with 
4+  2-f- 

Zr  and  Ca  ions  and  that  anion  vacancies  were  created  to  preserve 
electrical  neutrality.  The  present  experiment  was  performed  to  find 
direct  X-ray  evidence  for  the  proposed  model. 


The  intensity  of  X-ray  diffraction  lines  vas  measured  on 
a  Norelco  diffractometer  using  Cu  radiation.  The  composition  studied 
contained  l4  mole  per  cent  CaO  and  was  heated  at  1400  °C  for  a  week,  after 
sintering  at  2000 °C  for  2  hours.'  The  structure  factors,  were  then 

calculated  by  the  relationship: 


Flbi.thu)]  =  KJ 


(i) 


where  T  ,  is  the  observed  intensity,  p  is  the  Lorentz-polarization 
obs 

factor,  M  is  the  multiplicity  factor  and  K  is  a  scale  factor.  The 
observed  structure  factors  were  then  compared  with  the  calculated  values. 
These  calculations  were  based  on  two  models:  (l)  all  cations  randomly 
located  on  the  cation  lattice  sites  with  compensating  anion  vacancies, 
and  (2)  completely  filled  anion  sites  with  the  excess  cations  located 
in  l/2,  l/2,  l/2  interstitial  positions,  which  were  considered  to  be 
the  most  probable  sites.  The  calculated  structure  factors  were  obtained 
by: 

j  271  iCbu^+kv,,, (2) 

'co/c.chU)  =  C 

and 

(3> 

where  JL  ,  is  the  tabulated  i.onic  structure  factor;  B  is  a  temperature 
/fobs 

factor  and  u,  v  and  w  the  coordinates  of  the  ions.  The  temperature  factor 

-l6  2 

used  here  was  2.3  x  10  cm  . 

The  results  obtained  are  listed  in  Table  II.  In  general, 
the  anion  vacancy  model  fits  the  observed  values  better.  The  disagreement 
factor,  R,  was  computed  according  to: 


p  |  Fjjbs,  Fca/c. 

K  =  DF0U 

and  the  values  obtained  were: 


1.  R  =  6.2#  for  the  anion  vacancy  model. 

2.  R  =  11.2#  for  the  cation  interstitial  model. 
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These  results  favored  the  anion  vacancy  model  over  the  cation 
interstitial  model.  However,  the  existence  of  a  mixture  of  both  models 
could  not  be  completely  excluded  on  the  basis  of  these  results. 

3.  Electrical  Conductivity 

Electrical  conductivity  was  measured  at  1000  cps  on  all  specimens 
which  were  heated  at  l400°C  for  one  week  after  sintering  at  2000° C  for 
2  hours.  Measurements  were  made  on  a  General  Radio  impedance  bridge  Model 
650A  with  an  external  null  detector.  Platinum  electrodes  were  applied  in 
the  form  of  paste  to  the  cylindrical  specimens  which  were  then  heat  treat¬ 
ed  at'~i400 °C .  The  specimen  was  spring  loaded  between  two  platinum  blocks 
to  which  pla,tinum  lead  wires  were  attached.  The  whole  assembly  was  heated 
in  a  platinum  wound  tube  furnace  at  a  heating  rate  of  10 °C  per  minute  and 
measurements  were  made  at  20 °C  intervals  from  500  to  l400°C.  The  results 
are  presented  in  Figure  3  where  log  conductivity  is  plotted  versus  recip¬ 
rocal  absolute  temperature.  Between  and  about  1100 °C  the  conduc¬ 

tivity,  CT  (ohm  ^-cm  "*")  can  be  represented  by: 

<r  =  AexP  ^  (5) 

where  Q  is  the  activation  energy,  k  is  the  Boltzmann  constant,  T  is  the 
absolute  temperature  and  A  is  a  pre- exponential  term.  The  conductivity 
-temperature  data  for  the  various  cubic  solid  solutions  are  given  by  the 
following: 

for  13  mole  °jo  CaO 
for  15  mole  $  CaO 
for  16  mole  $  CaO 
for  18  mole  $  CaO 
for  20  mole  $  CaO 


1.26  x  103  exp  :1^ 


1.52  x  10 J  exp 


2.17  x  10J  exp 


1.71  x  10J  exp 


1.37  x  10J  exp 


-1.17 

kT 

-1.06 

kT 

-1.30 

kT 

-1-35 

■  kT 


ohm  "'‘-cm  ^ 


.  -1  -1 

ohm  -cm 


,  -1  -1 
ohm  -cm 


ohm  ^-cm  ^ 


.  -1  -1 
ohm  -cm 


H 


The  results  for  the  .15  mole  per  cent  CaO  composition  may  be  com- 
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pared  with  those  of  Hund  and  of  Kingery  et  al.  For  this  composition, 
Hund  obtained  an  activation  energy  of  1.21  e.  v.  and  Kingery' s  data  can 
be  expressed  by  (T  =  1.50  x  103  exp  (-1.26/kT).  Thus,  the  activation 
energy  for  conduction  obtained  by  different  investigators  is  in  fair 
agreement.  However,  the  conductivity  values  differ  widely.  Iiund's  speci¬ 
mens,  which  were  prepared  at  1500°C  and  therefore  were  perhaps  more  porous 
than  Kingery 's  and  our  specimens,  have  the  lower  conductivity.  Though 
Kingery's  specimens,  like  ours,  were  prepared  at  2000 °C,  his  conductivity 
values  are  lower  than  those  reported  here.  This  may  possibly  be  due  to 
contact  resistance  effects  as  a  result  of  employing  spring  loaded  platinum 
blocks,  whereas  painted  platinum  electrodes  were  used  in  the  present  study. 
The  conductivities  obtained  by  different  investigators  at  1000<’c  for 
ZrQ  CaQ  01  Q,-  are  compared  below: 


Conductivity,  ohm  1cm  1 

Investigator 

4  x  10-3 

Trornbe  and  Foex11 

2.7  x  10“3 

Volchenkova  and  Pal/guev1° 

2.2  x  10-3 

Hund® 

_2 

,  12 

2.3  x  10 

Kingery  et  al. 

2.6  x  10-2 
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Rhodes  and  Carter 

3-3  x  10"2 

Present  Work 

The  improved  conductivity  of  the  present  material  may  be  attribut¬ 
ed  to  several  factors;  higher  purity  materials,  higher  density  and  possibly 
better  electrodes. 

The  present  results  show  that  with  increasing  CaO  content,  the 
activation  energy  increases  and  the  conductivity  decreases  for  any  given 
temperature.  However,  the  conductivity  maximum  reported  by  Trornbe  and 
Foex11  and  by  Volchenkova  and  Pal'guev10  at  15  mole  per  cent  CaO  was  not 
observed  in  the  present  work.  The  variation  in  the  pre- experimental  term 
with  CaO  content  observed  in  the  present  study  is  not  considered  significant. 

IT 

The  electrical  conductivity,  (T ,  of  an  ionic  solid  can  be  written 


as : 


<T  =  YL  ini 


(6) 


where  rr  is  the  density  of  carriers  of  type  i  with  effective  charge  of 
q  and  mobility M  .  In  the  introduction,  evidence  has  been  cited  showing 
that  the  cation  and  electronic  contributions  to  the  total  conductivity 
of  the  fluorite  phase  in  the  ZrOg-CaO  system  are  very  small.  Therefore, 
anion  vacancies  dominate  the  conduction  mechanism.  Hence,  equation  (6) 
may  be  rewritten  as: 


{T  =  n  q M  (7) 

where  all  the  symbols  are  now  referred  to  oxygen  ions.  The  concentra- 
tion  of  anion  vacancies  is  fixed  by  composition,0  i.e.  each  molecule  of 
CaO  added  introduces  one  anion  vacancy.  It  may,  therefore,  be  expected, 
as  in  the  case  of  impure  alkali  and  silver  halides,  that  the  conductivity 
increases  essentially  linearly  with  defect  concentration.  However,  the 
present  data  (Figure  3)  show  that  the  conductivity  actually  decreases  as 
the  CaO  content  is  increased  from  13  to  20  mole  per  cent.  This  is  shown 
more  clearly  in  Figure  4,  where  the  logarithm  of  conductivity  was  found 
to  decrease  with  CaO  content  at  an  approximately  linear  rate,  the  rate 
of  this  decrease  decreasing  with  increasing  temperature.  This  behavior 
can  arise  from  either  or  both  of  the  following:  (1)  a  decrease  in  charge 
carrier  concentration  and  (2)  a  decrease  in  carrier  mobility,  with  in¬ 
creasing  CaO  content. 

In  the  solid  solutions  under  discussion,  the  defect  concentration 

(Ca  ions  and  anion  vacancies)  is  quite  high.  Under  these  conditions,  the 

charges  tend  to  balance  themselves  locally  so  that  the  probability  for  the 
2+ 

formation  of  Ca  ion  -  anion  vacancy  complexes,  or  probably  larger  clusters, 
is  quite  high.  Lidiara^  has  shown  that  at  low  concentrations  (C<  ^  1$) 
the  degree  of  association  increases  with  increasing  defect  concentration. 
However,  the  number  of  free  carriers  continues  to  increase  with  increasing 
defect  concentration  but  undoubtedly  at  lower  than  a  linear  rate.  Thus 
the  simple  association  theory  does  not  account  for  decreasing  conductivity 
through  reduced  carrier  concentration. 

The  following  model  is  proposed  to  account  for  the  observed  depen¬ 
dence  of  conductivity  on  CaO  content.  The  oxygen  ion,  which  is  the  charge 
carrier,  has  to  pass  between  two  metal  ions  to  reach  an  adjacent  anion  site. 


4+  4+  2+ 

These  metal  ions  may  be  two  Zr  ions,  one  Zr  and  one  Ca  ion,  or  two 

2+  2+ 

Ca  ions.  In  as  much  as  the  Ca  ion  (0.99  A)  is  approximately  25$ 

larger  than  the  Zr^h  ion  (0.78  A),  it  is  expected  that  the  energy  required 

•  2+  . 

for  an  oxygen  ion  to  pass  between  two  Ca  ions  will  be  largest  while  that 

A+ 

for  the  case  of  two  Zr  ions  is  the  smallest.  Each  oxygen  ion  is  surrounded 

by  four  metal  ions  in  the  fluorite  lattice.  As  the  CaO  content  increases 

2+ 

from  13  to  20  mole  per  cent,  the  probability  of  having  one  Ca  ion  as  a 

nearest  neighbor  to  an  oxygen  ion  increases  from  3^$  to  4l$  and  the  prob- 

2+ 

ability  of  having  two  Ca  ions  as  nearest  neighbors  to  an  oxygen  ion 
increases  from  7*7$  to  17$.  Therefore  the  activation  energy  for  conduc¬ 
tion  increases  and  consequently  the  conductivity  at  a  given  temperature 

decreases  with  increasing  CaO  content.  ■  The  above  calculation  is  based 

2+  2+ 
on  a  random  distribution  of  Ca  ions.  If  an  ordering  of  the  Ca  ions 

and  anion  vacancies  takes  place,  as  discussed  In  the  next  section,  the 

2+ 

probability  of  an  oxygen  ion  passing  a  Ca  ion  in  reaching  an  anion  vacancy 
increases  and  therefore  the  conductivity  is  lowered. 

Kingery^  obtained  good  agreement  between  the  diffusion  coefficient 
calculated  from  his  conductivity  data  using  the  Nernst-Einstein  relation 
and  the  measured  diffusion  coefficient  for  oxygen  ions  in  Zr^  85Ca0.15°1.85* 
This  fact  was  the  basis  for  attributing  the  electrical  conductivity  to 
oxygen  ion  mobility.  Anion  diffusion  coefficients  calculated  using  conduc¬ 
tivities  of  the  present  study  and  the  Nernst-Einstein  relation  are  higher 
(by  about  50$)  than  values  observed  by  Kingery.  Two  possible  factors,  if 

applicable,  would  tend  to  bring  the  calculated  and  observed  diffusion 

17 

coefficients  closer:  (l)  correlation  factor  which  for  fluorite  structure 
has  a  value  of  0.6-5*  and  (2)  defect  ordering  which  may  take  place  at 
lower  temperatures  during  the  diffusion  experiments  which  may  make  the 
use  of  lower  conductivities  of  the  ordered  state  more  appropriate’. 

In  the  temperature  range  1000  to  1200 °C,  the  log  O'  vs  l/T  curves 
show  a  change  in  slope,  the  magnitude  of  which  increases  with  increasing 
CaO  content  (Figure  3).  This  is  believed  to  be  a  true  behavior  of  the 


*  The  correlation  factor,  f,  is  defined  by:  ^  -  =  f  LkT/A/%?) 

where  !?.,  07  ,  and  q.  are  the  diffusion  coefficient,  electrical  con- 

i'  L  *  i 

ductivity  and  charge  of  i  ions,  respectively.  N  is  the  concentration 
of  lattice  sites  of  a  given  kind.  The  other  symbols  have  their  usual 
significance. 


material  and.  not  related  to  equipment  or  measurement  techniques.  The 
change  in  slope  arises  from  an  order- disorder  transition  as  discussed 
in  the  next  section. 

4,.  Order-Disorder  Transition 

Specimens,  sintered  at  2000 °C  for  2  hours  and  annealed  at  l400°C 
for  one  week,  shewed  only  the  cubic  fluorite  phase  by  X-ray  diffraction. 
However,  on  further  annealing  the  same  specimens  at  1000°C  for  one  week, 
extra  diffraction  lines  appeared  which  were  not  present  in  the  l400°C 

I 

annealed  samples.  The  extra  lines  could  not  be  accounted  for  by  monoclinic 
or  tetragonal  ZrOg,  or  CaZrO^ •  Thus  a  decomposition  of  the  cubic  phase  at 
1000  °C  can  be  ruled  out.  No  second  phase  was  observed  with  the  petrographic 
microscope. 

The  indexed  X-ray  patterns  of  the  20  mole  per  cent  CaO  sample 
heated  at  l400°C  and  at  1000*0  are  given  in  Table  III.  The  pattern  for 
the  1000  °G  sample  was  indexed  by  taking  a  cell  edge  four  times  the  basic 
cubic  cell  edge.  The  superstructure  lines  disappeared  when  the  specimens 
were  reheated  to  l400°C  indicating  a  reversible  order- disorder  transition. 

Electrical  conductivity  versus  temperature  data  for  a  20  mole  per 
cent  CaO  specimen  annealed  at  1000 °C  for  a  weak  are  shown  in  Figure  5*  At 
temperatures  below  1100 °C,  those  specimens  heated  at  1000  C  exhibited  a 
lower  conductivity  than  those  heated  at  1400  C.  Beyond  1100  C,  the  1000  C 
heat  treated  specimens  exhibit  a  change  in  slope  of  the  log  (f  versus  l/T 
curve  and  become  more  conductive.  On  reheating  the  same  specimen,  the 
log  (T  versus  l/T  plot  gives  a  straight  line  with  a  displacement  of  the 
curve  to  a  higher  conductivity  at  lower  temperatures  and  joining  the  first 
curve  at  higher  temperatures.  This  fact,  coupled  with  the  X-ray  diffraction 
data  were  the  evidence  of  an  order- disorder  transition  in  these  solid  solu¬ 
tions.  The  ordered  phase  was  less  conductive  o.ni  was  transformed  into  the 
disordered  phase  at  about  1100 °C.  The  disordered  phase  was  preserved  meta- 
stably  and  exhibits  a  higher  conductivity  on  reheating.  At  1000 °C,  the 
conductivity  of  the  disordered  phase  is  approximately  50$  larger  than 
that  of  the  ordered  phase  of  ZrQ  .  _  ^0  _  qq  for  this  particular  heat 

treatment. 


The  same  effect  was  observed  for  all  compositions;  however,  it 
was  more  pronounced  at  higher  CaO  concentrations  which  indicates  that  a 
higher  percentage  of  the  ordered  phase  is  formed  at  higher  CaO  concen¬ 
trations  . 

The  available  powder  X-ray  diffraction  data  do  not  permit  an 
understanding  of  the  nature  of  the  ordering  taking  place  in  these  solid 
solutions.  However,  the  fact  that  the  ordered  state  exhibits  lower  con¬ 
ductivity  than  the  disordered  state  must  imply  that  anion  vacancies  assume 

ordered  sites  instead  of  being  randomly  distributed.  The  following  is 

2+ 

therefore  proposed  as  a  plausible  ordering  mechanism.  The  Ca  ions  assume 

ordered  cation  sites.  This  is  facilitated  by  their  large  concentration  (13 

to  20  mole  per  cent)  and  the  low  cation  diffusion  rates'*^  (D  2+  ^  10 cm^  sec 
—  IT  P  _  I 

and  Dgr^  ^  1°"  cm  sec  in  ZrQ  g^  CaQ  ^  0^  g^  at  1000°C*).  The  ordered 
Ca^+  ions  attract  anion  vacancies  to  the  nearest  neighbor  positions  due  to 
coulombic  forces.  Thus,  the  anion  vacancies  become  ordered  leading  to 
decreased  ionic  conductivity.  Confirmation  of  this  model  may  be  sought 
via  neutron  diffraction  studies  of  ordered  solid  solutions.  The  order 
-disorder  transformation  observed  here  is  somewhat  analogous  to  the  transi¬ 
tion  from  a  disordered  fluorite  to  an  ordered  pyrochlore  structure  reported 
in  the  ZrCt,-  rare  earth  oxide  systems.^ 

IV.  SUMMARY 


1.  In  specimens  prepared  from  high  purity  ZrOg  and  CaCO^  at  2000°C, 
fluorite-type  solid  solutions  exist  from  between  12  and  13  mole  per  cent  to 
between  20  and  21  mole  per  cent  CaO- 

2.  X-ray  intensity  data  are  consistent  with  the  structural  model 

I4.4.  P-l- 

with  Zr  and  Ca  ions  completely  filling  the  cation  lattice  sites,  and 
oxygen  ion  sites,  equal  in  number  to  the  mole  per  cent  CaO  added,  being 
left  vacant. 

3.  In  the  temperature  range  5 00  to  1000°C,  the  electrical  conduc¬ 
tivity  can  be  represented  by  an  Arrhenius -type  expression  (F=  exp  (-Q/kT) 
where  the  activation  energy,  Q,  increases  with  increasing  CaO  content. 
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*  Obtained  by  extrapolating  Rhodes  and  Carter's  data. 
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4.  The  logarithm  of  conductivity  decreases  with  increasing  CaO 
content  essentially  at  a  linear  rate,  the  rate  of  decrease  decreasing 
with  increasing  temperature. 

5-  The  X-ray  diffraction  patterns  of  samples  annealed  at  1000°C 

exhibit  super- structure  lines  which  disappear  upon  heating  the  specimens 

to  1400 °C-  At  temperatures  below  about  1100 °C,  the  conductivity  of  these 

specimens  heated  at  1000 #C  is  approximately  50$  lower  than  those  heated 

at  1400 °C.  Above  1100 "c,  both  conductivity  curves  coincide.  This  behavior 

has  been  interpreted  in  terms  of  an  order- disorder  transition,  involving 
2+ 

the  Ca  ions  and  the  associated  anion  vacancies. 

ACKNOWIEDGEMENT 

The  authors  are  grateful  to  R.  R.  Heikes,  P.  H.  Sutter,  W.  G.  Carlson, 
R.  Ruka  and  J.  Bauerle  for  helpful  discussions. 


sr 


II 


REFERENCES 

1.  K.  Kiukkola  and  C.  Wagner,  J.  Electrochem.  Soc.  104,  379  (1957)* 

2.  J.  Weissbart  and  R.  J.  Ruka,  J.  Electrochem.  Soc.  lOg,  723  (1962). 

3.  E.  Ryshkewitch,  Oxide  Ceramics,  Academic  Press,  New  York,  i960, 

Chapt.  II.  5. 

4.  R.  S.  Roth,  J.  Am.  Ceram.  Soc.  39>  196  (1956). 

5.  J.  Stocker,  Bull.  Soc.  Chem.  France,  Jan.  1961,  p.  76. 

6.  R.  Collongues,  M.  P.  Jorba  and  J.  Lefevre,  Bull.  Soc.  Chem.,  France, 

Jan  1961,  p.  70. 

7.  P.  Duwez,  F.  Odell,  F.  H.  Brown,  Jr.,  J.  Am.  Ceram.  Soc.  ^5,  107  (1952). 

8.  F.  Hund,  Z.  Physik.  Chem.  199,  142  (1952). 

9.  A.  Dietzel  and  H.  Tober,  Ber.  deut  Keram.  Ges.  30,  47,  71  (1953)* 

10.  Z.  S.  Volchenkova  and  S.  F.  Palguev,  Trans.  Inst.  Electrochem.  No.  1 

Consultants  Bureau,  New  York,  1961,  p.  97* 

11.  F.  Trombe  and  M.  Foex,  Compt.  rend,  236,  1783  (1953)* 

12.  W.  D.  Kingery,  J.  Pappis,  M.  E.  Doty  and  D.  C.  Hill,  J.  Am.  Ceram. 

Soc.  42,  393  (1959). 

13.  W.  H.  Rhodes  and  R.  E.  Carter,  Presented  at  the  Am.  Ceram.  Soc.  meeting 
April,  1962,  Abstract  in  Bull.  Am.  Ceram.  Soc.  April  (1962). 

14.  J.  Weissbart  and  R.  J.  Ruka,  Electrochemical  Soc.  Extended  Abstracts 
No.  44,  Oct.  1961. 

15.  J.  B*  Nelson  and  D.  P.  Riley,  Proc.  Phys.  Soc,  (London)  57,  l60  (1945)* 

16.  See,  for  instance,  H.  Lipson  and  W.  Cochran,  The  Determination  of 
Crystal  Structure,  G.  Bell  &  Sons  Ltd.  (ixrndon)  1953* 

17.  A.  B.  Lidiard,  in  Handbuch  der  Physik  (S.  Flugge,  ed. )  Vol.  XX, 
Springer- Verlag,  Berlin,  1956. 


II 


TABLE  I 

SFECTROCHEMICAL  ANALYSIS  OF  ZrOg 


ELEMENTS  DETECTED 

ESTIMATED  CONCENTRATION 

Fe 

0.009  -  0.09$ 

Mg 

0.0003  -  0.003$ 

Cu 

0.0002  -  0.002$ 

Si 

0.0006  -  0.006$ 

TABLE  II 


X-RAY  DIFFRACTION  INTENSITY  DATA  FOR  POWDERED 
ZrO . 86  Ca0.i4  °i.86  SPECIMENS 


calc 


l  k  1 

f~dbs 

ANION 

VACANCY 

MODEL 

CATION 

INTERSTITIAL 

MODEL 

111  * 

137-74 

138.95 

126.64 

200 

78.49 

84.76 

89.83 

220 

127.96 

138.89 

147.20 

311 

94.71 

93.^5 

85-17 

222 

76.73 

66.28 

70.25 

400 

82. 4i 

93-09 

98.65 

331 

71.45 

67.47 

61.49 

420 

53-17 

50.38 

53-40 

422 

69.42 

68.02 

72.09 

511 

(333) 

53-01 

50.42 

45.96 

440 

50.36 

51.32 

54.39 

531 

40.89 

39-58 

36.07 

442 

(6oo) 

38.06 

29.95 

31.73 

620 

37-55 

39.80 

42.18 

533 

30.57 

30.16 

27-48 

TABLE  III 


X-RAY  DIFFRACTION  PATTERNS  FOR  Zr^^  CaQ  20  °l  80 

OBSERVED  d  VALUES  CALCULATED  d  VALUES 

Annealing  Temperature 

l400°C  1000°C  a  =  5.142  A  a  =  20. 568  A 


d,  A 

d,  A 

hkl 

a,  a 

hkl 

d,  A 

3-64 

440 

3.636 

3.28 

620 

3.250 

2.96 

2.97 

111 

2.963 

444 

2.968 

2.57 

2.57 

200 

2.571 

800 

2.571 

2.01 

10,2,0 

(862) 

2.017 

1.8l6 

1.815 

220 

1.8182 

880 

1.8182 

1.715 

12,0,0 

(884) 

1.7143 

1.549 

1.549 

311 

1.5503 

12,4,4 

1.5503 

1.483 

1.483 

222 

1.4845 

888 

1.4845 

1.284 

1.315 

4oo 

1.2855 

16,0,0 

1.2855 

1.179 

1-179 

331 

1.1796 

12,12,4 

1.1796 

1.150 

1.150 

420 

1.1498 

16,8,0 

1.1498 

1.050 

1.049 

422 

1.4095 

16,8,8 

1.4095 

0.9968 

0.9883 

511 

(333) 

0.9898 

20,4,4 

(12,12,12) 

0.9898 

0.9074 

0.9084 

440 

0.9091 

16,16,0 

O.9091 

0.8684 

0.8694 

531 

0.8690 

20,12,4 

O.8690 

0.8561 

0.8567 

600 

(442) 

0.8572 

24,0,0 

(16,16,8) 

0.8872 

O.8125 

0.8129 

620 

0.8130 

24,8,0 

0.8130 

0.7839 

0.7841 

533 

0.7842 

20,12,12 

0.7842 

(a)  12  mole  %  CaO,  (monoclinic 
Zr02  grains  are  shown  by  arrows) 


(b)  16  mole  %  CaO 


Fig.  1-Microstructure  of  Zr02  -  CaO  ceramics 
(2000  °C  -  2  hours,  1400  °C  -  one  week) 
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Fig.  4-Composition  dependence  of  conductivity  of  Zr02  -  CaO 
solid  solutions  at  indicated  temperatures 
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Fig.  5-Temperature  dependence  of  conductivity  of  Ca^  £q  Z^q  g qO^  gQ 
(2000°C  -  2  hrs.,  1400°C  -  one  week,  1000°C  -  one  week)  Showing 
order  -  disorder  transition 
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ORDER-DISORDER  AND  IONIC  CONDUCTIVITY  IN  CUBIC  ZrOg-CaO* 

E.  C.  Subbarso  and  P.  H.  Sutter 

In  the  fluorite-type  cubic  Zr^Ca^g  solid  solutions 

(0.12  <  x  <  0.22),  the  Zr^  and  Ca2+  ions  completely  fill  the  cation 

sublattice  while  a  fraction,  x,  of  anion  sublattice  sites  are  -vacant. 

(l)  (2-k) 

The  work  of  Wagner'  '  and  others'  '  has  established  that  the  electrical 
conductivity  of  these  solid  solutions  is  due  entirely  to  the  migration  of 
the  oxygen  ion  vacancies.  This  system  has  attracted  considerable  research 
effort  recently. ^  ^ 

Tien  and  Subbarao^  have  recently  found  an  electrical  con¬ 
ductivity  change  that  depends  on  the  heat  treatment  of  the  specimen 
and  have  interpreted  this  behavior  as  the  result  of  an  order-disorder 
transition.  They  also  found  superlattice  lines  in  the  x-ray  diffraction 
pattern  associated  with  the  ordered  state.  Additional  evidence 
concerning  the  order-disorder  phenomenon  and  the  associated  ionic  con¬ 
ductivity  changes  are  presented  here. 

A  cylindrical  specimen  of  which  118(1  been 

sintered  at  2000°C  and  painted  with  platinum  paste  electrodes,  was  used 

•Supported  by  ARPA  through  the  U.S.  Office  of  Naval  Research 
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for  the  electrical  conductivity  measurements  at  1000  cps.  The  upper 
line  in  the  log  a  vs  l/T  plot  (Fig.  1)  represents  the  conductivity  in 
the  disordered  state  as  measured  during  rapid  heating  (2°C  per  minute) 
of  the  sample  which  had  been  quenched  subsequent  to  a  long  anneal  at 
l400°C.  The  lower  line  was  measured  during  cooling  at  2/3°C  per  minute 
after  a  64  hour  anneal  at  1000°C  that  had  resulted  in  the  establishment 
of  the  ordered  state  with  its  lower  conductivity.  Cycling  between  the 
ordered  and  disordered  states  is  possible  merely  by  repeating  the 
appropriate  heat  treatment  so  that  the  effect  is  clearly  due  to 
recoverable  structural  changes. 

The  conductivity  changes  with  time  during  the  disorder-order 
transition  were  also  studied.  The  specimen  was  obtained  in  the  dis¬ 
ordered  state  by  annealing  at  l400°C  in  the  conductivity  measuring 
apparatus  which  was  then  rapidly  cooled  to  about  500°C  by  withdrawing  it 
from  the  hot  zone  of  the  furnace.  The  furnace  was  then  stabilized  at 
the  desired  annealing  temperature  and  the  sample  holder  was  reinserted 
fully  into  the  hot  zone.  This  resulted  in  a  rapid  rise  to  the  isothermal 
annealing  temperature,  at  which  the  conductivity  was  measured  as  a 
function  of  time.  The  results  at  885°C  are  shown  in  Fig.  2.  About 
half  of  the  conductivity  change  occurs  in  100'  minutes,  but  th"e 
rate  of  change  becomes  slower  as  the  transition  continues.  Hence,  not 
surprisingly,  the  decay  is  not  an  exponential  process  with  a  single 


relaxation  time 
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The  diffusion  coefficient  for  the  cations  in  this  system  is 
many  orders  of  magnitude  less  than  that  for  the  oxygen.'  '  The  measured 
transition  rate  seems  to  he  too  slow  for  the  process  to  he  oxygen  dif¬ 
fusion  controlled  since  the  diffusion  length  would  be  about  0.1  mm  for 
a  time  of  100  minutes.  However,  it  is  consistent  with  a  cation  dif¬ 
fusion  controlled  reaction  over  distances  of  a  few  lattice  spacings. 

This  indicates  that  the  conductivity  change  is  associated  with  cation 
migration  over  distances  so  small  that  appreciable  composition  in¬ 
homogeneities  cannot  be  established.  The  importance  of  cation  ordering 
is  also  shown  by  the  fact  that  the  observed  x-ray  superlattice  lines  are 
primarily  determined  by  the  cation  arrangement  due  to  their,  .much  higher 
scattering  power. 

According  to  the  present  results  and  those  of  Tien  and  Subbarao, 
the  conductivity  change  is  interpreted  as  an  order-disorder  phenomenon 
rather  than  the  precipitation  of  a  second  phase.  Order-disorder  requires 
only  a  rearrangement  of  the  atoms  by  means  of  a  few  atom  jumps  without 
any  composition  change.  Precipitation  implies  the  growth  of  regions  of 
a  second  phase  which  have  a  composition  different  from  the  high  tempera¬ 
ture  material,  while  regions  which  remain  in  the  fluorite  phase  would 
decrease  in  calcium  content.  Conductivity  changes  with  increasing  Ca 

content  in  the  ZrOg-CaO  system  have  recently  been  attributed  to  a 

(5) 

decomposition  of  the  cubic  phase,'  '  which  is  apparently  supported  by 
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}  (j  Q) 

earlier  x-ray  studies.  •  *  ■  Similar  arguments  may  be  advanced  to 
account  for  the  conductivity  changes  with  time  that  are  discussed  here. 
It  is,  however,  concluded  that  these  changes  arise  from  an  order- 
disorder  effect  since:  1}  According  to  the  phase  diagram, the 
likely  second  phases  are  monoclinic  (tetragonal)  ZrOg  and  CaZrO^. 

There  is  no  trace  at  all  of  these  phases  in  the  x-ray  patterns.  On 
the  contrary,  the  patterns  show  superlattice  lines  that  are  expected 
for  an  ordered  arrangement  of  the  lattice  into  a  larger  unit  cell. 

The  intensity  of  these  superlattice  lines  increases  with  increasing 
Ca  content  as  does  the  magnitude  of  the  conductivity  change.  2)  The 
kinetics  of  the  transformation  are  consistent  with  a  cation  migration 
of  only  a  few  atom  jumps  and  are  too  fast  to  allow  time  for  segrega¬ 
tion  into  phases  of  different  composition.  3)  The  lattice  parameter 
of  the  fluorite  lattice  does  not  show  any  change  of  composition  with 

S 

ordering  although  precipitation  of  a  second  phase  would  imply  the 
depletion  of  the  Ca  content  in  the  remaining  fluorite  phase.  4)  No 
phase  separation  was  observed  in  polished  sections  under  the  microscope, 
even  in  samples  which  exhibited  superlattice  diffraction  lines. 
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Conductivity  vs.  temperature  of  ZrQ  32  ^ao  i8°l  82 


Figure  1 


Isothermal  change  in  conductivity  at  885° C 
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CONDUCTION  PROPERTIES  OF  CdgNb^  AND  Fb^NbgOg^ 
R.  Mazelsky  and  W.  E.  Kramer 


Introduction 


Pyrochlore  is  a  mineral  of  the  approximate  composition  NaCaNbgOgF. 

Materials  of  this  structural  type  became  of  interest  when  ferroelectric 

( i  p) 

properties  were  observed  in  CdgNbgO^.. '  *  Cook  and  Jaffe  further  reported 
a  defect  pyrochlore  of  the  composition  Pb,  rNbo0^-  ...  Jona.  Shirane  and 

l.p  c  b.p 

(^*) 

Fepinsky  clarified  the  structure  of  CdgNbgO^.  They  reported  a  network 
of  NbOg  corner  sharing  octahedra  with  the  seventh  set  of  oxygens  and  the 
cadmium  ions  occupying  open  spaces  in  the  lattice.  These  relatively  "un¬ 
attached"  ions  lead  to  the  possibility  of  ionic  conductivity.  The  same 

situation  exists  for  Fb  cNb_Or  _  except  that  Borne  of  the  lead  and  oxygen 

J-O  O#? 

sites  are  vacant.  On  the  other  hand  electronic  conduction  is  possible  as 
a  result  of  the  variable  valency  exhibited  by  both  lead  and  cadmium  oxides. 
(CdO  via  Cd(m)  interstitials  and  FbO  with  its  positive  two  or  four  valence). 
An  attempt  will  be  made  to  elucidate  the  conduction  type  in  these  materials. 


Preparation 

The  materials  were  made  from  NbgO^  (99* 7$  minimum  purity)  CdO 
and  FbO  (Fisher  Certified).  For  the  purpose  of  comparison  cadmium  and 
lead  oxides  were  prepared  from  99 •999-*#  metal  by  means  of  solution, 
precipitation  and  careful  oxidation  to  the  monoxide. 
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The  oxides  were  weighed  to  yield  a  product  of  the  appropriate 
composition.  The  reactants  were  mixed  thoroughly,  compacted  into  a 
cylinder,  and  prefired  at  800°  for  24  hours.  Weight  losses  of  less  than 
1  part  per  thousand  were  observed.  The  pellet  was  placed  in  a  platinum 
cylinder  and  put  into  a  quartz  tube  which  was  previously  sealed  at  one 
end.  A  small  excess  of  CdO,  equivalent  to  the  weight  loss  observed,  was 
added.  The  quartz  tube  was  evacuated  and  back  filled  with  oxygen  to 
about  1/4  atmosphere.  The  quartz  tube  was  then  sealed.  The  tube  was 
heated  near  the  melting  point  of  CdgNbgO^  -  approximately  l4O0°  -  for 
3  hours.  A  similar  procedure  was  used  for  the  preparation  of  Fb^  j-NbgO^ 
A  final  firing  temperature  of  approximately  1290°  was  used. 

The  resulting  products  were  rounded  on  the  edges  -  indicative 
of  incipient  melting.  They  were  quite  crystalline  and  orange  in  color. 
Small  pieces  were  chipped  off  for  x-ray  identification  photographs  and 
symmetrical  cylinders  were  cut  out  for  measurements.  All  samples  used 
for  measurement  had  a  density  95-4$  of  theoretical. 

Experimental 

Using  the  procedure  described  above  a  solid  solution  series 
of  the  type  (l-x)  CdgNbgO^  -  x  Fb^  ^.NbgOg  ^  was  prepared.  Chemical 
analyses  were  made  on  samples  where  x  =  0.15  and  x  =  0.85  and  the 
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compositions  as  reported  agree  with  the  analysis  within  experimental 
error.  The  results  are  shown  on  Table  1. 


.85  CdgNbgO,-,  - 

Table  1 

•15 

.15  Cd2Nb207  - 

Fbl.^ 

Theor. 

Exp. 

Theor. 

Exp. 

Fb 

8.72 

8.39 

44.86 

43.4 

Cd 

35.76 

35-6 

5.73 

6.36 

Nb 

34.78 

35-5 

31.55 

32.0 

A  complete  range  of  solutions  appears  to  exist,  the  lattice  parameter 
varying  nonlinearly  with  composition  as  shown  on  Fig.  1.  The  powder 
diffraction  photographs  were  taken  with  a  Phillips  unit  having  114.6  mm 
diameter. 

Resistance  measurements  as  a  function  of  temperature  were  made 
using  an  impedance  bridge.  The  pellets  used  were  cylindrical  specimens, 
the  ends  of  which  were  plated  with  platinum  to  decrease  contact  resistance. 
Spring  loaded  platinum  electrodes  were  used..-  The  sample  and  holder  were 
heated  in  a  resistance  furnace  in  air.  The  results  are  shown  on  Fig.  2. 

The  activation  energies  calculated  from  the  In  p  vs  1000/T  are  shown  on 
Fig.  3.  A  minimum  in  the  activation  energy  is  observed  at  about  "x"  =  0.15* 
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If  one  plots  isothermal  resistivity  vs  composition,  at  all  the  tempera¬ 
tures  between  800°C  and  1150°C,  a  resistance  minimum  is  observed  at 
"x"  *  0.10  with  a  broad  maximum  over  the  center  range  of  composition 
(x  =  0.3  -  0.8),  the  resistivity  dropping  back  to  approximately  that  of 

Pb1.5Bb^,6.5  at  x  =  °‘9‘ 

Cylindrical  pellets  were  used  to  measure  EMF  with  a  Pg/p^  =  5 
where  P  and  P.,  were  pressures  of  oxygen  on  opposite  sides  of  the  pellet. 
Pig.  4  is  a  schematic  of  the  holde'r  used.  The  thermocouples  are  used  for 
temperature  measurement  on  both  sides  of  the  specimen.  The  platinum  legs 
on  both  sides  are  further  used  as  lead  wires  to  measure  EMF.  The  holder 
proper  is  made  of  lavite,  the  steel  rods  being  of  a  length  such  as  to 
keep  the  springs  outside  the  furnace.  The  specimen  faces  are  platinum 
coated  to  ensure  contact  between  the  thermocouple  and  the  sample.  The 
reaction  that  occurs  is  assumed  to  be  Og  +  4e  -*  2CT  and  the  percent 
oxygen  ion  conductivity  was  based  on  the  calculated  value.  Due  to 
the  high  resistance  of  the  sample,  the  sensitivity  of  the  EMF  measure¬ 
ments  was  low  below  850°C.  Measurements  were  taken  at  850°C  and  900°C. 
The  results  are  summarized  on  Table  2. 
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Table  2 

£  Ionic 

850°C 

Conduct. 

900°C 

Porosity  ($) 

Ca01>5  -  Zr02>85 

94 

~  7 

CdgNbgO^, 

57 

12 

~  6 

.95  CdgNb207  -  .05  Fb1>5m>2°6. 

5  U 

7 

~  8 

.9  Cd2Nb207  -  .1  Fb^NbgO^ 

10 

5 

~  5 

O  CdgNb 207  -  .5  Pbi.5Nb2°6.5 

0 

0 

~  5 

0 

0 

1 

Resistivity  measurements  as  a  function  of  oxygen  pressure 
were  made  on  CdgNbgO,^  and  Fb^  ..NbgOg  Oxygen  pressures  of  1. 5  mm 
to  760  mm  were  used  and  the  log  pressure  was  plotted  versus  log 
resistivity.  Cd^NtgO^  showed  little  if  any  dependence  on  oxygen 
pressure  below  900°C.  At  1030°C  a  variation  in  resistance  was  ob¬ 
served  with  pressure,  but  the  results  were  not  conclusive.  Measure¬ 
ments  of  Fb^  ,-NbgOg  j.  at  800°C,  900°C,  and  1000°C  are  shown  on 
Fig.  5*  The  resistivity  is  observed  to  vary  with  the  pressure  of 
oxygen  to  the  l/4  power. 

Discussion 

At  present  writing  final  conclusions  have  not  yet  been  made. 
However,  certain  observations  can  be  made. 
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1.  Below  1000°C  CdgNbgO^  seems  to  be  predominately  an  ionic 
conductor.  This  is  verified  by  both  EMF  measurements  and  the  independence 
of  resistivity  with  oxygen  pressure. 

2.  Fb.  cNbo0/-  c  is  an  electronic  conductor.  The  mechanism 

1*5  ^  De  P 

by  which  the  electronic:  carriers  are  provided  has  not  yet  been 
clarified.  However,  some  insight  can  be  gained  by  the  dependence  of 
resistivity  on  oxygen  pressure. 

3.  The  variation  of  resistivity  with  composition  is  not  clear. 
The  minimum  observed  at  about  .9  CdgNbgO^  -  .1  Fb^  ^Nb^Og  can  be  ex¬ 
plained  by  the  presence  of  electronic  conduction,  but  the  broad  maximum 
observed  over  the  center  range  is  not  clear  —  particularly  since  EMF 
measurements  show  these  compositions  are  electronic  conductors.  The 
dependence  of  resistance  on  oxygen  pressure  has  not  yet  been  determined 
for  these  compositions,  but  this  work  will  be  done. 
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Captions 

-  Lattice  Parameter  vs 

"x" 

(mole  fraction  Pb^  ^NbgOg 

-  P/PQ  vs  1000 /T 

Identification  No. 

Pc 

> 

"x"  (mole  fraction  Pb1  ^Nbg0g 

1 

4.4  x 

io-6 

0.00 

.4 

2 

1.2  x 

10 

0.05 

3 

1.8  x 

10-- 

0.15 

4 

3.9  x 

10-5 

0.50 

5 

5-5  X 

10"5 

0.70 

6 

3.1  x 

10”5 

o.Qy 

7 

3.4  X 

10"  5 

1.00 

-  Activation  Bnergy  vs 

"x" 

(mole  fraction  Pb^  j-Nb^Og  ^ ) 

Fig.  4  -  EMF  Measurement  Holder 

A.  Sample 

B.  Gold  0  rings 

C.  Pt  -  Pt  10?t  RL  Thermocouple 

D.  Quartz  delivery  tube 

E.  Ceramic  Rod 

F.  Steel  Rod 

G.  Spring 

Fig.  5  -  Log  resistivity  vs  log  of  the  oxygen  pressure  at  the  temperatures. 
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Abstract:  The  transport  properties  of  mixed  valence 
semiconductors  are  studied.  In  addition  their  dependence  on  the 
magnetic  properties  is  considered.  It  is  assumed  that  the  conduction 
process  is  a  thermally  activated  one  in  which  the  charge  carriers 
jump  from  site  to  site.  The  activation  energy,  A  H« ,  required  for 
a  jump  is  calculated  assuming  that  a  jump  can  only  take  place  when  ' 
the  atomic  configuration  surrounding  a  charge  carrier  and  that 
surrounding  a  neighboring  site  are  energetically  equivalent.  This 
activation  energy  is  related  to  the  force  constants  and  in  turn  to 
the  differences  In  size  of  the  cation  in  the  two  valence  states.  This 
trapping  mechanism  is  shown  to  be  more  plausible  than  the  long  range 
polarization  effects  previously  proposed.  A  calculation  of  the  Seebeck 
coefficient,  Ot  ,  on  this  model  shows  that  the  contribution  due  to  the 
activated  process  is  normally  small  and  of  the  order  of  —  AHf  . 

20  eT 

Magnetic  effects  give  rise  to  a  small  contribution  to  cc  and  slightly 
modify  the  temperature  dependence  of  the  electrical  conductivity.  No 
discontinuities  in  the  transport  properties  are  expected  at  the  Curie 
temperature.  In  order  to  understand  the  properties  at  high  carrier 
concentrations  it  is  found  necessary  to  modify  the  concept  of  free 
carriers,  and  in  addition  consideration  must  be  given  to  determine  the 
number  of  sites  accessible  to  the  carrier. 
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Resume :  On  a  etudie  les  proprietes  de  transport  des  semiconducteurs 
de  valence  melangee.  De  plus  leur  dependance  des  proprietes  magnetique 
est  consideree.  On  suppose  que 'la  conduction  prend  place  peer  une  diffusion 
thermique  activee  par  laquelle  le  porteur  saute  d’un  atome  a  1' autre.- 
L’energie  d' activation  pour  un  saut,  AHC  ,est  calculee  en  supposant  qu*un 
saut  ne  peut  prendre  place  que  si  la  configuration  atomique  entourant  le 
porteur  et  celle  d'un  site  voisin  sont  energiquement  equivalantes .  Cette 
energie  d' activation  est  liee  aux  coefficients  locaux  d'elasticite  et  a 
son  tour  a  la  difference  du  taille  du  cation  dans  les  deux  etats  de  valence. 

On  montre  que  ce  mecanisme  de  pi^gement  est  plus  plausible  que  celui  de  la 
polarization  a  grande  distance  propose  anterieureraent.  Un  calctil  du 
pouvoir  thermoelectrique,  oC  ,  "base  sur  ce  modele  montre  que  la  contribution 
a  OC  due  au  processus  active  est  normalement  petite  et  de  l'ordre  de 
—  ^  Hc„  .  Les  effets  magnetiques  occasionnent  une  petite  contribution 

*o  eT 

a  OC  et  modifient  legerement  la  dependance  de  la  temperature  de  la 
conductivite  electrique.  Aucune  discontinuite  des  proprietes  de  transport 
ne  s' attend  a  la  temperature  de  Curie.  Afin  de  comprendre  les  effets  a  de 
hautes  concentrations  de  porteurs  il  faut  modifier  la  conception  de 
porteurs  libres  et  en  oivtre  determiner  soigneusement  le  nombre  de  sites 
actuellement  accessibles  aux  porteurs . 
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Une  Etude  des  Proprietes  de  Transport  des 
Semiconducteurs  de  Valence  Melang£e  (*) 

R.  R.  He  ike#1 
A.  A.  Maradudin 
■  and 

R.  C.  Miller 

INTRODUCTION 

Dans  les  semiconducteurs  non-ferromagnetiques  de  valence 

melangee  le  mecanisme  du  transport  de  la  charge  paraft  etre  une 

diffusion  thermique  activee  par  laquelle  les  porteurs  sautent  d'ua 

atome  a  1' autre. Un  cas  typique  est  NiO.  Dans  ce  cas  on  peut 

introduire  des  porteurs  en  remplaipaut  une  petite  quantite  de  Ni  pa;r 

Li.  On  suppose  alors  que  les  conditions  de  la  neutrality  de  charge 

occasionnent  la  formation  de  Ni  .  La  charge  supplementaire  sur  lea 
+3 

Ni  ,  le  trou,  est  done  le  porteur.  Sous  certaines  conditions  il 
est  possible  que  les  porteurs  soient  pieges  par  la  distortion  creee 
par  eux-mene .  Or  le  trou  ne  peut  se  mouvoir  que  si  l'on  a  un 
enlevement  approprie  de  cette  distortion.  On  a  suppose  que  la 
frequence  des  sauts  peut  s'ecrire  dans  un  tel  cas^"^ : 

^"R.  R.  Heikes  and  W.  D.  Johnston,  J.  Chem.  Phys.  2 6,  582-7  (1957) • 

*Cette  oeuvre  a  ete  effectuee  pendant  que  cet  auteur  a  se  joume  quelques 
mois  au  Centre  National  de  la  Recherche  Scientifique  a  Bellevue. 

(*)  Ce  travail  a  ete  subventionne  en  partie  par  Office  of  Naval  Research, 


Denartment  of  the  Wav'/ 
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-A  Gc/kT, 

>>  =  v0  z  e  ,  (i) 

ou  est  l'energie  libre  qui  est  necessaire  pour  enlever,  de  la  . 
maniere  precisee  ci-dessous,  la  distortion  entourant  le  trou  en 
question.  Normalement  on  trouve  que  la  frequence  a  une  grandeur 

tres  proche  des  modes  optiques.  7L  est  le  nombre  d'atomes  voisins 
auquel  le  trou  peut  sauter. 

Auparavant^  on  a  calcule  les  proprietes  de  transport  en 
utilisant  (l).  On  a  montre  que  la.  conductivite  electrique  depend 
essentiellement  de  */T,  tandis  que  le  pouvoir  thermoelectrique,  OL  , 
comprend  trois  termes. 


ou  c  est  la  concentration  des  porteurs  libres  et  Ce  la  concentration 
des  sites  accessibles  aux  porteurs.  (On  verra  ci-dessous  qu'il  faut 
modifier  la  signification  normale  de  "porteurs  libres"  a  des  concentrations 
>Z7..)  Le  jjremier  terme  depend  fortement  de "la  dynamique  du  mouvement 
du  porteur.  Le  deuxieme  terme  est  simplement  le  changement  d ' entropie 
quand  un  porteur  supplementaire  est  ajoute  et  piege.  Le  troisieme  terme, 
1' entropie  de  melange,  ne  depend  que  des  concentrations  C  et  Co  • 

2j(.  p.  He ikes  and  R.  W.  Ure,  Jr.,  Editors,  Thermoelectricity:  Science 
and  Engineering,  Interscience  Publishers,  New  York,  19 6l. 
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Dans  cet  article,  nous  voulons  developper  1' etude  anterieure 
dans  deux  directions.  D' abord,  on  etudiera  l'effet  des  proprietes 
magnetiques  du  compose  sur  les  propri6t4s  de  transport.  En  particulier 
on  examiners  ce  qui  se  passe  au  voisinage  de  la  temperature  de  Curie,  Tt 
Deuxiemement  on  etudiera  le  mecanisme  du  mouvement  d'un  porteur  afin  de 
mieux  comprendre  les  quantites  AGC  ,  A  ,  et  AS  qui  se  trouvent  dans 
(1)  et  (2). 

Enfin  on  se  servira  de  ces  calculs  pour  construire  un  modele 
assez  complet  des  proprietes  de  transport  des  semiconducteurs  de  valence 
melangee . 

II.  La  Conductivity  Electrique 

En  utilisant  (l)  on  neglige  le  fait  que  la  probabilite  d'un 

saut  pourrait  dependre  de  1 ' orientation  relative  des  spins  des  deux 

atomes  en  question.  En  fait  Anderson  et  Hasegawa'J'  ont  mont.r£  que  la 

probabilite  en  depend  assez  fortement.  Dans  la  limite  classique  de 

tres  grand  spin,  ils  montrent  que  la  probabilite  d’une  transition  est 

proportioxmelle  a  cos®  ou  0  est  1' angle  entre  les  deux  spins  en 

2 

question.  Ici  meme  si  la  distorsion  est  enlevee,  il  se  peut  que  la 
transition  ne  prenne  pas  place.  On  doit  done  remplacer  (l)jpar 


^P.  W.  Anderson  and  A.  Hasegawa,  Phys.  Rev.  100,  675  (1955)* 
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-AGt/lcT 


\ 

ou 


\>  *  Vo  zfm  e 


f  (T)  =  Z  n^./w  . 

0 


(3) 


(4) 


n#  est  le  nombre  de  paires  d'atomes  vole ins  dont  les  spins  font  un 
angle  0  j  "Pg  est  la  probabilite  d'une  transition  dans  ce  cas  et  w 
est  le  nombre  global  de  paires  d'  atomes .  Pour  simplifier  le  probleme 
nous  supposons  qu'un  atome  donne  ne  peut  avoir  que  deux  etats  de  spin, 
soit  parable le,  soit  antiparallele  a  une  direction  donnee.  Alors 
nous  prenons  P0  ■  I  ,  P(  a  0.(Quoique  cette  hypothese  puisse  changer 
la  Pome  exacte  du  resultat,  elle  ne  changera  pas  les  conclusions.) 

Etudions  maintenant  la  conductivity  llectrique  dans  un 
compose  antiferromagnetique  qui  a  un  reseau  cubique  simple.  Prenons 
la  disposition  des  spins  ou  le  reseau  cubique  est  decompose  dans  deux 
reseaux  face  centres,  chacun  ferromagnetique  mais  se  trouvant  anti¬ 
parallele  l'un  a  l'autre.  Nous  nous  servons  de  l'equation  d'Einstein. 

Nee1  D 


c r  = 


JcT 


(5) 


ou  c  est  la  concentration  de  porteurs  libres,  N  le  nombre  d’atomes  par 
unite  de  volume,  et  D  la  constante  de  diffusion  des  porteurs.  ^ 
s '  ecrit 


D  =  saw 


(6) 
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ou  £  est  un  coefficient  geometrique  (e-gal  a  l/6  dans  ce  cas),  a  est 
la  longueur  d'un  saut  et  V*  est  donne  par  (4).  Nous  supposons  que  le 
porteur  ne  peut  sauter  qu'aux  premiers  voisins.  II  ne  reste  plus  qu'a 
determiner  la  forme  de  la  f onction  f ( T) .  Dans  notre  reseau  cubique 
simple,  la  disposition  des  spins  est  telle  que  chaque  paire  de  premiers 
voisins  se  trouvent  antiparallele  a  T  —  0 .  (Nous  n'allons  pas  considerer 
le  cas  ferromagnetique  car  il  parait  que  les  porteurs  ne  se  piegent  pas 
dans  ce  cas-la.)  Si  l'on  considers  le  transport  de  la  charge  entre 
deux  plans  (100)  dans  une  direction  [100],  il  est  evident  que  chaque 
saut  change  le  sous-reseau  du  porteur.  Ainsi,  la  forme  de  la  fonction 
•f  peut  s'ecrire 

f  *  <*.  +■  Fx  •  (7) 

Ot,  et  jS ,  representent  la  prohahilite  qu'on  trouve  un  atome  sur  le 
sous-reseau  I  avec  le  spin  par allele  ou  antiparallele  respectivement  et 
une  direction  donnee.  L‘ indice  2  refere  au  deuxieme  sous-reseau. 

Puisque  le  sous-reseau  2  est  antiparallele  au  sous-reseau  I  ,  nous 
voyons  que  0^  a  |3(  et  jj^s  oc(.  Done,  nous  pouvons  refere  tout  au 
sous-reseau  I  et  recrire  (7)  comme 

f  =  2a,  =  2«(i-cc)  (8) 

ou  nous  avons  supprime  1' indice  I  .  Si  nous  exprimons  ■(*  en  fonction 
de  M/M0  )  ou  Mo  est  1' aimantation  d'un  sous-reseau  a  Ts  0  et  M  est 
l'aimantation  a  une  temperature  T  ,  il  est  clair  que 


Z  oc  -  I  • 


Done 


f-  if  -  (an- 

De  1' equation  (5)  nous  pouvons  ecrire  la  conductivity  comme 

-4G«At 


\  ( 

m  \M  e 

do) 

(J~  —  constante 

•  “  l 

■ 

hJ  j 

T 

Nous  pourrions  proceder  en  utilisant 

une  fonction  de  Brillouin  pour 

M  ;  pourtant,  la  siraplicite  de  notre  modele  ne 

le  justifie  pas. 

Supposons  plutot  que 

(Mf  »  ,  - 

/  T  \ 

X 

1 

(11) 

IM./ 

1  V 

\  -p  / 

ou  1 t  est  la  temperature  de  Curie. 

Ceci  est  une  approximation 

raisonnable  d'une  fonction  de  Brillouin.  Nous  avons  done 


-AGcAT 


(12) 


<r  -  T  e  au-dessous 

,  -A  Ge/k,T  -r 

ty  ~  _L  Q  au-dessus  |c  , 

T 

II  est  evident  que  cette  dependance  de  la  temperature  sera  cachee  si 
AG  c  est  suffisament  grand. 
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III.  Le  Pou voir  Thermoelectrique 

II  nous  faut  maintenant  examiner  le  me'canisme  du  transport  de 
la  charge.  A j ou'tons  un  porteur  supplementaire  a  un  atome  du  reseau 
tout  en  tenant  fixe  tous  les  autres  atomes.  (Dans  la  suite  nous 
apellerons  souvent  un  atome  avec  une  charge  supplementaire,  un  hote 
charge^,  et  les  autres  atomes,  les  hotes.)  Puis  si  nous  liberons  les 
atomes  hotes,  ils  se  deforment  a  cause  de  cette  charge  supplementaire. 

Par  consequent,  l'energie  litre  du  systeme  decro£tra.  Nous  denommons 

par  AG  le  changement  de  l'energie  litre  quand  une  charge  supplementaire 

H 

est  ajoutee  et  pie'gee  a  un  atome  donne'.  Afin  que  la  charge  puisse  Be 
mouvoir,  nous  supposons  que  les  configurations  atomiques  entourant  la 
position  de  la  charge  supplementaire  (position  initiale)  et  la  position 
a  laquelle  elle  saute  (position  finale)  dolvent  etre  ldentiques.  Or, 
les  positions  pourraient  devenir  identiques  en  enlevant  completement  la 
distortion  a  la  position  initiale.  D'une  mainiere  semblahle  on  pourrait 

^  N  ^  \ 

creer  une  distortion  a  la  position  finale  egale  a  celle  de  la  position 
initiale.  Ou  plus  generalement,  on  pourrait  enlever  partlellement  la 
distortion  a  la  position  initiale  et  creer  une.  distortion  semblable  a 
la  position  finale. 

Or,  si  le  travail,  depense  a  la  position  initiale  pour  mouvoir 
une  charge,  est  dc-sirne  par  W ,  snivant  la  reference  1,  nous  ecrivons 
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ou  -  ~  (SAG^/^T ).  AS^^est  le  changement  de  la  contribution 

magnetique  a  l'entropie  due  a  1' addition  d'une  charge  supplementaire. 
Dans  notre  modele  simplifie  ce  terme  est  identiquement  zero,  car  nous 
n'avons  permis  que  deux  etats  de  spin,  ne  dependant  pas  de  la  presence 
ou  de  1' absence  d'une  charge  supplementaire.  Toutefois,  dans  un  cas 
reel,  1* addition  d'une  charge  supplementaire  accroit  le  nombre  d' etats 
de  cet  atome,  done,  l'entropie.  Mais  a  T*  0  ,  ou  nous  avons  un  etat 
antiferromagnetique  parfait,  AS  j^q*0,  PuiS(lue  1' addition  d'une 
charge  supplementaire  n ' accroJt  pas  le  nombre  de  configurations. 
Cependant,  au-dessus  de  Tt  ,  ou  tous  les  etats  de  spin  sont  egalement 
probable,  il  est  aise  de  montrer  que 


AS 


MAC 


=  Cvt 


2Sa,+  1 

2  S ,  +  I 


ou  s  represente  le  spin,  1' indice  I  se  refere  aux  hotes  charges  et  2 


aux  hotes. 

Pour  un  compose  tel  que  Li^Ni^^O,  Sj  »  3/2,  S|  s  '  /  AS  MAG  “ 
Cette  contribution  au  pouvoir  thermoelectrique  au-dessus  de  sera 
done  de  l'ordre  de  ZSyv/'Q.  On  s' attend  a  ce  que  cette  contribution 
soit  zero  a  T"0  et  croisse  jusqu'  a- <w"25ja»/*C  a  ,  restant  constante 
au-dessus  de  7^  .  Ainsi,  cette  contribution  n'est  pas  importante  dans 


0.2?lc. 


m*. 
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V 


la  plupart  des  cas.  II  faut  remarquer  qu'on  ne  prevoit  pas  une 
discontinuity  d  . 

II  nous  faut  maintenant  etudier  la  contribution  provenant 
du  piegement  et  du  mecanisme  du  mouvement.  De  (13)  nous  voyons  qu'il 
faut  determiner  as*  et  W  -  Supposons  que  les  voisins,  entourant  la 
charge  supplementaire  se  deplacent  une  distance  V  dans  l'£tat  piege. 

Nous  introduisons  un  parametre  X,  qui  peut  varier  de  zero  a  I  .  Or, 
afin  de  laisser  le  porteur  se  mouvoir  il  faut  que  la  position  initiale 
et  la  position  finale  soient  identiques.  Nous  faisons  cela  en  enlevant 
partiellement  la  distortion  d  la  position  initiale  (c'est  a  dire  nouB 
laissons  les  voisins  se  deplacer  une  distance  XV  au  lieu  de  V)  et  en 
creant  la  meme  distortion  a  la  position  finale  (nous  deplacons  les  voisins 
entourant  la  position  finale  une  distance  XV).  Ainsi  le  travail  fait  pour 
laisser  mouvoir  le  porteur  est 


AG 


(HO 


F  /  A  / 

.  est  l'energie  litre  quand  les  atomes  entourant  l'hote  charge  ont  un 
/  r  *■ 

deplacement  XV.  est  l'energie  litre  quand  les  atomes  entourant  un 
hote  ont  un  deplacement  XV.  La  premiere  parenthese  represente  le 
changement  de  l'energie  litre  a  la  position  initiale,  tandis  que  la 
deuxieme  represente  le  changement  de  l'energie  litre  a  la  position  finale. 


V 
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p  s  constante  +  \  2.8  —— —  ( l  -  —  —  p"( _£.  \  fi  _  i.  £. 

rk  2  r°  i:  w  <*>••/  hv\^i  ('  ?$« 


Maintenant  recrivons  (3-5)  et  (3.6) 

X.  X4 

p,  s  constante  +  (  Y  ~  A  )  A  +  T  I) 

(Remarquons  que  AG0  *  F.s  constante  -  1  A  +TD) 

K.  l  Z 


(16) 

(17) 


~  constante 


+  ^  c . 


(18) 


k  * 

A,  C,  et  se  definissent  en  comparant  (3.7)  et  (l8)  avec  (15)  et  (l6). 
Done  on  peut  e valuer  AG  en  substituant  (17)  et  (l8)  dans  (l^). 

X1  \  >  a  ,  A  „ 


AG  •=  ( A  .  X)A  *  A  * 


(19) 


Nous  supposons  que  le  mecanisme  actuel  de  transport  soit  celul  qui 
requiere  la  valeur  minimale  d.e  AG.  II  est  aise  de  voir  que  cette 
condition  s'obtient  quand 

.  A 

k  =  - — r  ■  (20) 


A+C 


PuiE 


A  G. 


=  A 


c/A 


I  +  C/A 


(21) 


Nous  allons  identifier  cette  guantite  avec  I'energie  litre  d' activation 
pour  un  saut.  Puisque  AQC  ne  depend  pas  de  la  temperature  dans  notre 
modele  simplifie,  I'energie  d,' activation  pour  un  saut,  AHc,  est  egale 


a  A(j r . 
c 

Afin  de  calculer  w,  il  faut  avoir  le  changement  de  I'energie 
litre  a  la  position  initiale  seulement.  Evidemment  nous  avons 


w 


(22) 


car  on  depense  du  travail  en  deplacant  les  atomes  de  V  a  XV  tandis  qu'on 
recupere  du  travail  quand  les  atomes  se  deplacent  de  XV  a  zero  en 
l'atsence  d'un  porteur.  Or,  en  utilisant  (22)  et.les  equations 
precedentes,  nous  pouvons  evaluer  (13) 


e  ot 


c 

A 

($■•) 

IT  ( 

1*0 

I2 

(23) 


II  nous  reste  done  a  trouver  une  approximation  raisonnatle  pour  (23)  et 
(21).  Pour  faire  cela  nous  nous  servons  des  grandeurs  suivantes  pour 
les  parametres : 

(a)  On  montre  dans  1' appendix  qu'on  peut  ecrire 

=  M  cop  /H  , 


.  ,  / 
ou  00  e  est  une  frequence  moyene  efficace  du  reseau. 
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(b)  La  separation  d'equilibre  ^  dans  le  cristal  pur  est 
la  separation  a  laquelle  4s  =  0  •  Ainsi  4>'(r, »)  =  o, 

(c)  Nous  pouvons  developper  %(r)  en  serie  autour  de  , 

la  separation  a  laquelle  X'(rh  0.  Ainsi  X‘(r)  «  X'(r.)  +  (<■-  r,)X"(r,) 

D'ou 

p‘(r6)  =  (r0  -  r;  )%"(0  =  -Ar  X  (r, ) . 

>  (24) 

(d)  Puisque  nous  negligeons  les  contributions  anharmoniques, 

X"  M  =  XVr.)  =  co,z 


p"  s  £1  (  W*  ~  CO*  )  . 

1  M 

Avec  ces  donnees  nous  trouvons  pour  les  coefficients 

A  =  -  1) 

C  =  A  {  l  -  A(r-,)(i-  | )}  (£ 

D  =  3  Ic  {  R  -  I  -  2£R} 

ou  £  a  AV/r0  et  i  ^  .  On  voit  aue  le  parametre  £  est  une 
mesure  du  changement  du  rayon  a'un  hote  charge. 
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IV.  Discussion 

IIous  pouvons  maintenant  presenter  un  modele  assez  coraplet  qui 
nous  permet  de  comprendre  les  proprietes  de  transport  des  seraiconducteurs 
de  valence  melangee .  Pour'  mieux  definir  .notre  probleme  nous  donnons 
quelques  resultats  obtenus  sur  des  composes  typiques .  Dans  les  figures 

S  s  s  s 

1  et  2  nous  presentons  .le  pouvoir  thermoelectrique  et.  la  reslstivite 

electrique  pour  plusieurs  composes  du  systems  Li^Mn^  S.  Les  memes  donnees 

pour  le  systeme  Li  Mn.^  ^.Se,  se  trouve  dans  les  figures  3  et  it.  Enfin  pour 

qu'on  puisse  faire  une  comparaison  nous  montrons  les  resultats  pour  trois 

composes  diff'erents  de  meme  concentration  d'impuretes  (Figures  5  6). 

Dans  tous  -ces  examples  on  trouve  la  meme  structure,  cclle  du  sel  gemme ; 

toutefois  nous  considerons  que  ce  comportement  est  representatif  des  semi- 

ccnducteurs  de  valence  melangee.  On  pourrait,  par  exemple,  comparer  ces 

(it) 

donnees  avec  celles  d'un  perovskite. 

Afin  de  faciliter  la  discussion,  nous  idealisons  ces  donnees  dans 
la  figure  J.  Nous  pretendons  que  ces  courbes  montrent  toutes  les 
caracteristiques  necessaires . 

Examinons  maintenant  les  differentes  regions  de  ces  courbes. 

La  decroissance  exponentielle  de  la  reslstivite  avec  la  temperature  peut 
avoir  deux  origines:  (l)  la  liberation  des  porteurs  des  centres  accepteurs 

j|R.  C.  Miller,  R.  R.  He  ikes,  R.  Mazelsky,  J.  Appl.  Phy.  32,  2202  (1961). 
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(ou  donneurs)  avec  la  temperature.  Ceci  nous  donne  une  dependance  de  la 

_  q 

temperature  de  6  ou  Q  est  1' energie  de  liaison  entre  le  porteur  et 

le  centre.  (2)  La  dependence  exponentielle  avec  la  temperature  de  la 
frequence  d'un  saut.  La  decroissance  du  pouvoir  thermoelectrique  dans  la 

meme  region  de  temperature  provient  des  memes  origines.  II  est  evident 

If  Q 

que  la  liberation  de  porteurs  contribue  un  teme  de  la  forme  —  C:  ; 

.  eikT 

toutefois,  la  relation  entre  a  et  1' energie  d' activation  pour  un  saut 
n'est  pas  si  evidente.  En  \itilisant  (21)  nous  pouvons  recrire  (23)  comme 


e 

lc 


a  =  -  D  +  lilil  +  In  itii-  ■ 

Ic  kT  (§*,) 


(26) 


Nous  raisonnons  ci-dessous  que  C/A  doit  etre  tou jours  plus  grand  que  1' unite 
Ainsi  nous  voyons  que  ce  processus  donne  aussi  unterme  qul  decroit  avec  la 
temperature.  II  faut  remarquer  que  lorsqqe  C/A  est  egal  a  1' unite,  cet 
effet  n'existe  pas. 

A  des  temperatures  plus  hautes,  on  volt  que  la  resistivite  tend  a 
se  nivelcr.  Cela  peut  provenir  de  la  dependance  explicite  de  la  temperature 
trouvee  dans  (12).  Une  deuxi.eme  raison  pourrait  etre  l'epuisement  des 
porteurs  des  centres  donneurs  (ou  accepteurs).  Une  etude  de  (12)  sugg^re 
qu'en  traqarit  Tp/r,  au-dessous  de  Tc  et  If  /T  au-dess'us  on  trouverait  une 
ligne  droite.  Sur  la  figure  5  on  voit  le  resultat  pour  Li  .05M.95°-  Bien 
qu' on  n'obtienne  pas  de  resultats  tellement  excellents  dans  les  systemes  MnS 
el  I-lnGe  on  trouve,  en  effet,  que  la  courburc  tend  toujours  a  disparaitre. 
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le!  nous  pouvons  disposer  du  sujet  de  la  relation  entre  l'energie 
d' activation  pour  la  conduction  AH  et  Tt  .  (Remarquons  que  l'energie 
d' activation  pour  la  conduction  est  egale  a  AHC  +  Q.)  Sur  les  figures 

nous  voyons  que  certainement  AH  ne  subit  pas  un  changeraent  brusque  el  T  . 

c* 

C'est  a  cela  que  nous  devrions  nous  attendre  avec  notre  module,  car,  en 
derniere  analyse,  nous  exprimons  A  Hc  en  fonction  des  potentiels  d' interaction 
entre  les  atomes  qui  a  leur  tour  peuvent  s' exprimer  en  fonction  des 
coefficients  d'elasticite.  II  est  Men  connu  que  ceux-ci  ne  suMssent  pas 
de  changement  brusque  a  7^  .  (Les  changemen'ts  que  1'on  observe  sont 
attribues  au  mouvement  des  parois  de  Bloch.)  De  plus,  on  ne  s' attend  pas 
a  ce  que  Q.  subisse  un  changement  brusque  a  T^,  .  Cependant  du  fait  qu'il 
y  a  des  resultats  sur  quelques  composes  qui  montrent  un  changement  abrupt 
de  AH,  nous  devons  examiner  ce  point.  Wright  et  Andrews^ ^  et  Van  Houten^^ 
trouvent  une  discontinuite  dans  AH  au  voisinage  de  pour  NiO .  Aussi 
Heikes  et  Johnson  indiquent  une  telle  discontinuite  dans  CuO  aussi  bien 
que  WiO.  Wous  croyons  que  ces  discontinuites  ont  ete  mal  interpretees . 

D'abord  la  temperature  a  laquelle  la  discontinuite  se  produit  varle  de 

o  -  (6) 

plus  de  100  C  pour  plusieurs  composes  "purs"  de  WiO.'  'II  est  tres 

improbable  que  77  puisse-  varier  tant  dans  ce  cas.  De  plus,  cette 

discontinuite  disparait  quand  on  ajoute  une  petite  quantite  de  Li,  laquelle 

quaabitt  ce  change  pas  serieusement  le  comportement  magnetique.  Wous 

j 

^R.  W .  ;;ri  bt,  J.  P.  Andrews,  Proc.  Phys.  Soc.  A62,  446  (1949) . 

'  0.  Van  Won  ten,  J. .  Chon.  Phys.  Solids  17,  7  (i960). 
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attrlbuons  done  cet  effet  a  une  impurete  ou  plus  probablement  a  une 
propriete  de  la  lacune  de  Ni  trouvee  dans  KiO.  En  tout  cas  une  liaison 
avec  Tc  est  tres  invraisemblable . 

Quant  au  pouvoir  the  mo  e  le  c  t  r  i  que  a  des  temperatures  plus  hautes, 
nous  voyons  sur  la  figure  7  que  pour  des  concentrations  >  2*7o  n  commence 
a  augmenter.  Ensuite  nous  nous  attendons  &  ce  qu.'  il  se  nivele .  Par 
contre,  a  de  tres  basses  concentrations  ou  bien  cet  accroissement  n'a 
pas  lieu  ou  bien  il  est  tres  faible.  Comme  le  premier  terme  de  (26)  est 
independant  de  la  temperature  et  le  deuxieme  decroit  continuellement,  il 
faut  etudier  le  troisieme  pour  comprendre  cette  region.  Considerons 
d'abord  la  situation  quand  la  concentration  est  tres  petite.  Ici,  c  e8t 
s implement  la  concentration  de  porteurs  libres  determinfee  4  partir  du 
traitement  statistique  normal  et  c0  est  la  concentration  de  sites 
accessibles  aux  porteurs,  e'est  a  dire,  tous  les  sites  du  cristal  sauf  les 
impuretes.  (A  proprement  parler  on  devrait  soustraire  un  certain  nombre 
d'atomes  entourant  chaque  impurete;  pourtaat,  a  de  basses  concentrations 
ce  nombre  est  petit  devant  le  nombre  global.)  Ainsi  dans  cette  limite, 
le  troisieme  terme  decroit  continuellement.  Or,  a  des  concentrations 
plus  hautes  une  situation  entierement  differente  prevaut.  On  a  discute 
ce  cas  auparavant On  a  signale  que  le  porteur  prefere  se  trouver  au 
vois inage  du  nombre  maximum  d' impuretes  (Li  dans  not re  cas)  ce  que 
correspond  au  cas  de  l’energie  minimum.  Nous  voulons  poursuivre  ce  point 
de  vue  ici.  A  une  concentration  donnee  d' impuretes  il  y  aura  un  certain 
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pourcentage  des  cations  qui  n'a  aucune  impurete  comme  premier  voisin. 

II  y  aura  d'autres  qui  n'en  auront  qu'une,  d'autres  qui  n'en  auront  que 
deux,  etc.  Si  l'on  s-jppose  que  les  impuretes  se  distribuent  a.u  hs,zard 
partout  dans  le  reseau,  on  peut  calculer  ces  pourcentages  facilement. 

i 

On  trouve  les  resultats  pour  quelques  composes  du  systems  Li  Hi.  0  dans 
le  Tableau  I.  Or,  on  peut  montrer  que  pour  des  concentrations  au-dessus 
de  I'ordre  de  2$,  un  porteur  peut  se  mouvoir  partout  dans  le  cristal 
tout  en  restart  premier  voision  a:  une  impurete.  (lei  nous  ne  consid4rons 
que  les  sauts  a!  travers  un  oxygene,  c'est  a  dire  aux  deuxi^me  voisins. 

•Les  resultats  magnetiques  sur  ces  composes  nous  menent  a  supposer  que 
cette  transition  est  dominate.)  De  plus,  pour  des  concentration  au-dessus 
de  I'ordre  de  6fj>  les  porteurs  peuvent  trouver  des  chemins  o)i  ils  restent 
premier  voisin  a  deux  impuretes.  Voyons,  rnaintenout,  ce  que  veut  dire 
"porteur  libre"  dans  ce  cas.  Examinons  les  chiffres  pour  x  .03  sur  le 
Tableau  I.  Comme  il  vient  d'etre  dit,  les  cations  qui  ont  une  impurete 
comme  premier  voisin  constituent  des  chemins  continus  a  travers  le  reseau; 
cependant  les  cations  qui  ont  deux  impuretes  comme  premiers  voisins  n'ont 
pas  cette  propriete.  Done,  un  porteur  sur  un.  tel  site  est  piege  et  ne 
contribue  pas  a  la  conductivite.  Alors,  on  peut  envisiger  le  processus 
de  conductivite  de  la  maniere  suivante.  Les  trous  sont  excites  a  partir 
des  pieges  (les  sites  qui  ont  deux  ou  plus  de  deux  impuretes  comme 
premiers  voisins)  jusqu'au  greupe  d'etats  possedant  une  impurete  comme 
premier  voisin  ou  plus  haut,  jusqu'au  groupe  d'etats  ne  possedant  aucune 
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impurete  comme  premier  voisin.  ,  Ainsi  la  conductivite  a  lieu  sur  deux 
groupes  d'etats.  (Gertes  il  y  a  de  la  dispersion  entre  les  deux  groupes 
d'etats.)  De  plus,  nous  savons  que  le  pouvoir  themoelectrique  dans  ce 
cas  peut  s'lcrire. 


Oi 


<3~c  CX0 

+  °- 


(27) 


ou  crt  et  sont  les  conductivity  eleetriques  des  porteurs  sur  les 
sites  possedant  zero  ou  une  impure t!  respectivement  comme  premier  voisin. 

On  defini  (X0  et  Of,  dans  une  maniere  semblable .  Examinons  maintenant 
1' allure  d'  (X.  pour  x  =  .03-  A  de  basses  temperatures  o"0  est 
essentiellement  nulle;  OC  est  done  egal  a  Of,.  Or  & ,  decroit  avec  la 
temperature  a  cause  de  l'accroissement  des  porteurs  sur  les  sites 
possedant  une  impurete  comme  premier  voisin.  Quand  la  temperature  s'accrolt 
les  sites  de  plus  haute  energie  (n'ayant  aucune  impurete  comme  premier 
voisin)  deviennent  accessibles.  Alors  Oi  croit  et  enfin,  quand  toutes 
les  sites  sont  accessibles,  il  se  nivele. 

II  est  aise  de  montrer  que  1' analyse  donnee  ci  -  dessus  est 
correcte;  toutefois  on  trouve  de  grandes  difficultes  en  essayant  de 
poursuivre  le  calcul  quantitativement .  La  difficult!  actuelle  resulte 
du  fait  que  les  sites  possedant,  un  nonib  re  donne  d'impuretes  comme  premier 
voisin  ne  sont  pas  distribues  au  hazard,  plutot  ils  tendent  a  se  grouper. 
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Bien  que  nous  ne  discuterons  pas  la  region  de  concentration 
interraediare,  il  est  assez  facile  d'envisager  qualitativement  la 
transition. 

Ayant  discute  le  cadre  genlral,  examinons  la  signification 
du  calcul  de  la  Section  III  (les  details  se  trouvent  dans  1' appendix.) 
Nous  etudierons  ensuite  les  donnees  tenant  conrpte  du  modele  ci-dessus . 

Les  coefficients  de  (25)  qui  entrent  dans  (21)  et  (23)  sont  fonctions 
de  deux  pararaetres  £  et  "R  .  Pour  aider  notre  discussion,  il  est  peut- 
etre  plus  utile  de  nous  servir  de  C/A  et  el  ,  le  deplacement  d'un  atonie 
entourant  un  hote  charge,  comme  parametres  independents.  En  utilisant  la 
formulation  trouvee  dans  1*  appendix  on  pent  montrer  que 


2  £  eR 

ra  it 

De  (l'7)  et  (l8)  on  voit  que  C/A  est  le  rapport  de  l'energie 
requise  pour  cr4er  la  distortion  autour  d’un  hote  a  l'energie  requise 
pour  enlever  la  distortion  autour  d'un  hote  charge.  On  s' attend  5  ce 
que  dl/r0  soit  de  1'ordre  de  quelques  pourcents  -&t  que  C/A  3oit  pres  de 
l’unite.  Ainsi  nous  pouvons  recrire  (21). 


AG  —  A  , £/ ,  r*.  ^ .  =  constante  »  d  .  (29) 

C  ~  2  It  C/A  H 

Nous  recrivons  aussi  (23) 
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[Vu  (27)  il  est  clair  que  cette  equation  est  gross ierement  approcheej 
toutefois  elle  suffira  pour  la  suite,] 

Nous  voudrions  maintenent  nous  servir  de  (29)  et  (30)  pour 
etudier  des  composes  reels.  Deux:  objections  se  presentent  imraSdiatement : 

(1)  Nous  ne  considerons  que  d.es  interactions  entre  premiers  voisins, 

(2)  Notre  crista!  n’est  que  monoatomique .  Examinons  la  premier  objection. 

Nous  pretendons  que  le  mecanisme  de  piegement  dans  ces  composes 
n'est  actuellement  pas  la  polarization  a  grande  distance  propose  par 
Landau  mais  que  c'est.  un  effet  local  dependant  de  la  distortion  provenant 
pour  la  plupart  be  la  difference  du  rayon  de  l'ion  dans  les  deux  £tats  de 
valence.  Afin  de  supporter  cette  hypothilse,  nous  pourrions  mentiormer 
deux  fails  exp^rir.ientaux.  D'abord,  on  observe  ce  piegement  mo  me  ri  des 
concentrations  tres  hautes  d' impure  tes  (jusqu'a  25$).  De  plus,  la 
profondeur  du  piege,  determinee  par  l'lnergie  d' activation,  est  sensiblement 
constante  au-dessus  de  quelques  pourcents.  Certes,  a  une  concentration  de 
25$  il  n'est  pas  raisonnable  de  croire  que  le  piegement  est  du  ti  une 
polarization  k  grande  distance.  S'il  en  etait  ainsi,  nous  nous  attendrions 
a  trouver  un  changeraent  continu  de  l'energie  d'. activation  avec  la 
concentration.  Un  deuxieme  group  d' experiences  qui  porte  sur  cette 
question  est  1' elude  de  Li(ili,  Co)o.  ^  Ici,  les  effets  singuliers  ne  se 
laissent  s implement  comprendre  que  si  1'on  traite  separement  ehaque  saut 
possible,  c'est  a  dire,  qu'il  faut  considerer  que  le  porteur  peut  sauter 

C.  Miller,  2. 


Ti.  L.eibes,  W.  D.  Johnston,  J.  Chem.  Piiys.  Jl,  11 6  (1959) « 
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d'tm  Co  jusqu'a  un  Ni,  d'un  Wi  jusqu'a  un  autre  li,  etc.  Si  le  pi^gement 
etait  vraiment  un  effet  de  grande  distance  on  ne  eompterait  pas  distinguer 
de  tels  effets  locaux. 

Quoique  le  calcul  ait  ete  fait  pour  un  cristal  monoatomique 
nous  croyons  que  (29)  et  (30)  ecrites  en  fonctions  de  d  et  C//\  ont  une 
validite  plus  generale. 

II  faut  remarquer  que  nous  avons  traite  le  cristal  dans  le  cas 
de  1' approximation  harmonique.  (Nous  ne  trouvons  done  aucune  entrople 
d* activation. )  Nous  avons  neglige  les  effets  anharmonique  non  suelement 
a  cause  de  la  difficulte  du  traitement  mais  aussi  a  cause  du  fait  qu’il 
y  aurait  deux  parametres  de  plus,  diff  idles  a  a  valuer .  De  plus,  la 
simplicite  de  notre  module  ne  justifie  pas  tel  detail. 

II.  faut  maintenant  ‘brievement  examiner  les  facteurs  qui  reglent 
les  grandeurs  de  d  et  C//\  dans  tin  cristal  reel.  Dans  un  reseau  ionique 
nous  pretendons  que  1' effet  dominant  pour  la  determination  de  d  est  le 
rapport  du  rayon  de  1"  anion,  ,  a  celui  du  cation,  R(l  .  Nous  pouvons 
examiner  cette  question  sur  une  "base  purement  geometrique.  'Lorsque  le 

rapport  Vr,  est  suffisament  grand  pour  que  les  anions  soient  sur  le 

w  c. 

point  de  se  toucher,  il  est  evident  que  3„e  remplacement  du  cation  par  un 
autre  plus  petit  ne  peut  pas  aboutir  a  un  grand  deplacement.  Par  centre 

quandR.  /d  est  plus  petit,  1' interaction  anion-anion  n'empeche  pas  le 

A  & 

displacement  et  nous  nous  attendons  ainsi  a  un  plus  grand  deplacement. 
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Or,  nous  voyons  qua  AG  t  varie  comine  d  .  Ainsi,  on  corapte  que  l'energie 
d' activation  croit  de  MO  jusqu'a  MnO  car  le  rayon  du  cation  croit.  Aussi 
dans  la  serie  MnO,  MnS,  MnSe,  AGC  devrait  d£croitre  pour  la  mime  raison. 
Tout  cela  est  en  ton  accord  avec  1' experience.  Enfin,  voyons  si  des 
grandeurs  raisonnables  de  d  nous  donneront  le  propre  ordre  de  grandeur 
de  AGC.  Considlrons,  par  exemple,  Li  ^Wi  or0,  ou  AGt  ~  3ooo  cal. 

En  utilisant  (29),  (28),  (25)  et  prennant  CO„  ~  lO^sec"1,  f0  ~  H  *  1 0  cm, 
Ms  $  *  10  (9X1,  on  trouve  facilement  que  «•  ~  0.03.  Ceci  parait  Itre  un 

deplacement  raisonnable. 

Ensuite  discutons  la  grandeur  de  C/A  .  Afin  d'y  arriver  il 
faut  estimer  la  grandeur  de  E  (voir  (25)).  Certes,  quand  l'hQte  charge 
est  plus  petit  que  1'hote  meme,  on  s ' attendrait  et  ce  que  R  soit  moins  que 
1' unite  mais  du  meme  ordre  de  grandeur,  disons  entre  .5  et.l,  Ceci  veut 
dire  que  t/A  varie  de  1’ unite  jusqu'a  1.1.  Actuellement  on  peut 
conjecturer  que  C/f\  sera  toujours  plus  grand  que  1' unite  etant  independent 
du  rayon  de  1'hote  charge  car  Zener^a  montre  theoriquement  que  les 
coefficients  d'elasticite  decroissent  toujours  quand  on  ajoute  une  impurete. 
En  utilisant  les  grandeurs  ci  ~  dessus  pour  d  et  C/A  on  voit  que  D  [voir 

(25)]  doit  etre  moins  de  ~  .  En  genei'al  0x1  peut  done  negliger  cette 

l 

contribution  a  oi  ■ 

Einalcsnent  il  faut  faire  quelques  remarques  sur  l'energie  me  sure  e 
d* activation  pour  la  conduction.  Comme  nous  avons  deja  dit,  il  peut  y  avoir 


n 

,;C.  Zcnor,  Acta  Cryst. ,  (1949). 
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deux  contributions:  (l)  AGc  >  l'energie  d'activation  d'un  saut  et  (2) 
ou  Q.  est  l'energie  de  liaison  entre  le  porteur  et  1' impure t£ .  Maintenant 
nous  voulous  montrer  qu'  en  general  les  deux  effets  contribuent  mais  que  le 
premier  est  dominant  a  de  hautes  concentrations.  Considlrons,  &  nouveau, 


Li  „,.Ni  , r0 .  On  voit  sur  les  figures  5  et  6  que  l'energie  mesur&e 

.0?  .yy 

d' activation, AM  est  egale  a  3600  cal  et  que  la  contribution  a  Ot  dependante 


On  doit  done  supposer 


OH 
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Appendix:  Calcul  de  1' Energie  Libre 

Dans  le  calcul  suivant,  nous  traitons  un  reseau  cubique 
Bravais  qui  contient  N  mailles.  Nous  mettons  un  hote  charge  a  l'origine 
du  systems  des  coordonnees  et  nous  imposons  les  conditions  aux  limites 
periodiques  aux  deplacements  atomiques . 

De  plus  nous  supposons  que  les  atonies  hotes  interagissent 
deux  si  deux  avec  un  potentiel  .  D'une  maniere  semblable  nous 

supposons  que  l'hote  charge  interagit  avec  les  atomes  h&tes  ayant  un 
potentiel  Xlr).  L' energie  potentielle  du  cristal,  contenant  un  hi&te' 
charge  peut  s'ecrire 

V  =  *  A  4> 

ou 

<§>  -  i  £'  4>  (lrJ- 

1  tv 

a#  =  H'{X( (33) 

Dans  (32)  et  (33)  est  la  position  instantanee  du  ^eme  atome.  Le 

prime  sur  la  premiere  somme  exclut  les  termes  avec  e  -t  .  tandis  que 
le  prime  sur  la  deuxieme  somme  exclut  le  terme  a  0  .  Maintenant 

exprimons  sous  la  forme 


! 

1 
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r*  =  *lt)  +  (34) 

ou  *(0  est  la  position  vectorielle  du  {,:ieme  atorae  dans  le  cristal 
parfait  et  se  donne  par  la  formule 


?£(0  —  "^i  Q-i  +  +  2:3  '  (35) 

AV  ***•  ^ 

^  ^3  sont  des  noin:bre8  positifs,  negatifs,  ou  nuls  et 

Cl  i  CX,  Cl  »  s°nt  trols  vecteurs  non  coplanalres  caracteristiques  du 
cristal.  11(0  dans  ( 34)  represente  un  deplacement  du  atome  de 

sa  position  d'e quill! re  dans  le  cristal  parfait. 

Si  nous  developpons  l'energie  potentielle  en  puissances  de  U ({  ) 


nous  oltenons 


4=.  *  X  L  Ut')  M  a  f  «')+••• 


A<$>  =  A<i>  +  H'p«X)(u.,(0-  u.,o>i)  +  (3T) 

oc 

+  1  21'  u«(oi)(u.f(l)-u.,(o>) 
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Dana  (36)et  (37)  U*  U)  est  la  (X  -composante  du  deplacement  W.(0  . 
Dans  ( 37)  nous  avons  introduit  la  definition 

p  (<■)  =  Xlr)-  4>(r).  (38) 

Les  seuls  coefficients  non  nuls  dans  (36)  sont 

<#>.  =  i  1 I'  4>  0  Xtt)-  (39.) 

■U* 

=  (39M 

=  -  {  **  u)  xfU)  D^u)  +  S*f  V4>fr)} 

-  r'4vo  (39c) 

et, 

A4>.  =  H'  f  ('£'<>')  (to.) 

■e 
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P«U)  = 


-  2) 


ffr) 

i  r's  trt) 


=  xji)  DfCO 


(40b) 


*„  tt)  xf  (0  + 

(4oc) 


En  ecrivant  ces  re suit at s  nous  nous  servons  de  la  notation  suivante 


4>u)  ®  4? (u ^)i ) 

D4>(r)  =  -L  4>'  (r) 


(41) 


D*4>fr)  =  i  4>"  (r)  -  4=’Cr-^  . 

Dans  le  traitement  suivant  il  sera  convenable  de  pouvoir 
traiter  les  deplacements  U  (o)  et  U(£  )  (6fe)dans  (37)  sur  une  tase  egale. 
Nous  faisons  cela  afin  qu’il  ne  faille  pas  distlnguer  entre  les  derivees 
par  rapport  a  U<(o)  et  par  rapport  a  U*  (£  )  [it  0) . 

Nous  rlcrivons  done  (37)  formellement  sous  la  forme 

A  4?  =  A#.  +  Z  f,(e)u..tt)  * 

«  (te) 

Z  r 

le  coefficient  l.(w>  est  symmetrique  dans  l1  e change  de  liY 


ou 


avec 


* 

(  ^  ) .  Une  comparaison  de  (37)  et  (4-2)  nous  donne  les  relations  entre  les 
f  -coefficients  et  les  p  -coefficients.  Nous  voyons  que 


(i  )  •  =  p*  (t )  it  o 


(43a) 


*P«  (o) 


n't* uy  • 
•& 
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Si  nous  utilisons  (40b)et  (43b),  nous  trouvons  pour  notre  modele 

■fjo)  =  o  . 


(43c) 


Aussi,  nous  obtenons  la  relation  suivante  entre  (•(<')  et  '■ 

V<r)  =  ('■  P'fin- 

*"  m 

-  0- *«■)[*«.  p«,M'>  +  Sf. 

Si  nous  remarquons  que  (l )  =  p^,*  (£)  nous  voyons  que  les 

coeffecients  (iV)  donnes  par  (44),  satisfont  aux  relations 

•P«,  MM  =  f.«  MM)  = 

f  f  (45) 

Z  f«,  MM  =  O 


comme  il  le  faut. 


Maintenant  nous  ecrivons  U  (O  comme  la  somme  d'un.  deplacement 
statique  v(t)  et  un  deplacement  dynaanique  I id)  .  L’ energie  potentielle 


peut  s ' ecrire 


#  =  +  i  21  + 

AH'  r 

+  Z  #„/«')  frfMJVfM")  +... 

♦  f  21  r«M)  *,(*•)  +•■• 


(46a) 


A  #  =  A  +  r  21^ (*«*) 4-. 
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*  Jl  f-JO  SK(t)  +./L  $.U)v,(v)  +••• 


+  x  H  0«')  KM  s,w)  *■■  ■  ■  c**) 

-f*'  p  r 

*f>  - 

La  condition  de  l’equilihre  pour  le  cristal  est 

^77  ^  +  =  0  ■  £ $*,<“')  vr'm 

__  „  P  0*7) 

+  f««)  +  21  f„  Uv) *■■■ 

A'  '  r 

? 

Les  deplaeements  {.  V<8'  ^  ^  )  ,  c’est  a  dire,  les 

solutions  du  groupe  d' equations  (47),  decrivent  la  distortion  statiqu.e 
du  cristal  due  a  un  h8te  charge. 

Si  nous  introduisons  la  matrice  u  l1)  qui  est  1' inverse 

*r 

de  la  matrice  (<£')  c'est  a  dire 

*  y 

Nous  pouvons  alors  ecrire  (47)  dans  une  forrae  equivalence 

vj'YO  =  -  21  -  2ZG„ttu')ff1(tT)vrMu")* 

>  f  t  (lt9) 

En  resolvent  (4$))  par  iteration  nous  oh  tenons 

=  "  ZZ.  (W  -f  a  ^  )  + 

■C'  r  (50) 

P 


d  »  f  J*  :  ‘  1 
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Une  expression  explicite  pour  Gxp(H')  Ians  un  reseau  arbitraire 

V-  afjcJ)e,(U) 


Bi-avais  a  ete  trouve^et  elle  est 


Vu,)  * 


N 


nx(kj) 


(51) 


Dans  cette  expression  k  est  un  vecteur  dont  les  grandeurs  sont  distributes 
uniformement  dans  la  premiere  zone  Brillouin  du  cristal.  n  (U)  eat  une 
des  frequences  des  modes  no.rma.ux  pour  le  cristal  parfait  dan3  le  cas 
ou  la  masse  atomique  est  egale  a  1* unite.  L* indice  j  designs  la  polarization 
des  ondes  du  reseau  qui  correspondent  au  vecteur  d'onde  k  ,  tandis  que 
e(fcJ)  est  le  vecteur  unitaire  de  polarization  de  I'onde  (k  j)  . 

Uous  pouvons  maintenant  calculer  I’energie  statique  du  cristal 
deforme.  De  (46)  nous  voyons  que 

E,uti.  -  +  i  £4>f(U')v™tov?U')+- 

if  (52) 

+  Z  +  2  Z  Vw,)  <"ll)  +  ■ '  • 

Ot  /  *  $ 

au  second  ordre  des  deplacements  (  y^  U)}  .  Si  nous  nous  se:rvons  de 
la  condition  d’equilibre,  (4'f),  nous  trouvons  que 


Z.  {  4>Kf  UV)  +  hf  (w)}  v“’U )  v"W) 

=  -  H  tj*)v?u). 


41 


(53) 


En  utilisant  (52)  et  (53)  nous  ob tenons  enfin 


^P.  A,  Flinn,  A.  A.  Karadudin,  Annals  of  Physics  18,  8l  (1962). 
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E  stvh.c  '  4^  +*  +  i  v^u).  (54). 

* 

La  part  de  l'energie  litre  &’ Helmholtz  dependant  de  la 
temperature  provient  de  la  contribution  vibratoire.  Af in  ■  d*  ^valuer 
c.elle~ci  il  faut  determiner  les  frequences  des  modes  normaux  du  cristal 
deforme. 


Nous  commenfons  par  ecrire  1* equation  de  mouvement  du  cristal. 

M  (<§>  +  A  <$>) 

(»> 

-■  - 

Si  nous  prenons  une  solution  period! que  pour  ces  equations 


rj*’f-e)  e 


U  CO 


) 


(56) 


ou  •vo  est  independent  du  temps,  alors  les  equations  pour  les 

{STtt) )  sont 


o  y.  (o) 


u)  =  H  {  ^ (u‘]  <•  ^ 

r  M  m 


cm 

-  3E  u{l)  . 

V  r  p 

P 

Ainsi  les  carres  des  frequences  des  modes  normaux  du  cristal  deforme 
sont  les  valeurs  propres  de  la  matrice  ^o<a  • 
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/  »  3,  y 

Designons  les  valeurs  propres  &e  cette  matrice  par  j  00  j  \ 
On  pent  alors  ecrire  la  contribution  vibratoire  a  l*energie  libre  de 
Helmholtz  du  cristal  deforme  sous  la  forme 


3N 


=  kT  L  Ml } 


2  k  T 


(58) 


Dans  la  limite  de  hautes  temperatures,  qui  nous  interesse 
exclusivement  par  la  suite,  cette  expression  se  reduit  a 


!c  T 


3N 

L  * 

j 


/VC 


too; 

“kT 


(59) 


devaluation  exacte  de  (59)  pour  un  cristal  deforme  est 
extremement  difficile  et  on  n'y  est  pas  encore  arrive.  Done  nous  nous 
servons  d’une  methods  approchee  pour  evaluer  cette  expression.  Uous 
recrivons  (59)  sous  la  forme 


“  3  N  kT  .'U 


fa  co 

It 


(6o) 


Une  comparalson  de  (59)  et  (6o)  sert  a  definir  la  frequence  moyenne  CO 
Cependant,  pour  notre  but,  nous  pouvons  rapprocher  la  frequence  par' 


CO 


<M  ZSI 


T” , 


(6l) 
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\ 

ou 


GO1 


(62) 


Nous  faisons  ce  rapprochement  pour  deux  raisons.  D'abord,  nous  verrons 
que  la  quantite  CO  ^  est  facile  S,  evaluer.  Deuxidmement ,  ce  rapprochement 
est  exact  si  nous  utilisons  le  rapprochement  de  Debye  pour  le  spectre  des 
frequences. 

Maintenant  si  nous  nous  servons  du  theordme  que  la  somme  des 

i  6s  dig  s  \ 

n  '  puissances  des  valeurs  propres  d’une  matrice  est  egale  a  la  trace 

ifilTlG  / 

de  la  n  puissance  de  la  matrice,  nous  obtenons  immediatement 


CO 1  = 


3  N  L 


£  in) 


(63) 


_  J__  2Z  l  &■«(<<)+  h- (u)  }- 

i ,  l  DC  »<  ’  ■  Pig,-  J 

3NM  ; 

En  simplifiant  la  somme  dans  (64)  nous  trouvons  que 

H  --  ZI'4>«1W) 

4  *■  i> 


(64) 


IN 


X 

0< 


X 

oC 
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D'une  maniere  semblable,  avec  l'aide  de  (44)  nous  obtenons 


l 

O 4. 


=  2  H' {  p'Vxd))  +  p'(*«>)]  • 

-e  x(0  '  j 


(66) 


Afin  d'evaluer  les  expressions  que  nous  avons  obtenues, 
nous  aliens  etudier  un  reseau  cubique  face  centree.  De  plus  nous 
allons  supposer  que  les  potentials  4^>  et  X  (r)  n' interagissent 
qu' entre  premiers  voisins.  Ainsi  =  o  ,  ou  ro  est  la  separation 

equilibre  entre  premiers  voisins. 

y  /  * 

Nous  pouvons  evaluer  l'energie  de  la  distortion  du  reseau 
statique,  donnee  par  (54),  de  la  maniere  suivante.  Elle  est  donnee  par 

E.  =  1  Z^Ov-,'°VO 

z  JL  •  ' 


ip'  (67) 

r-  CO  Cl) 

‘  Es  +  Es  . 


*En  utilisant  (40)  et  (43)  pour  ■f«,u )  nous  pouvons  recrire 

(1) 


E  sous  la  forme 
s 


e‘"  =  -  x  (f^) 


(68) 


OU 


les  sornmes.  sur  (  et  ('  s’evaluent  pour -les  douze  premiers  voisins  de 


l'hote  charge,  qui  se  trouve  a  I'origine  des  coordonnees.  Si  nous  utilisons 

(q) 

le  tableau  des  fonctions  de  Green  de  Flinn  et  Maradudin,  nous  pouvons 
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e valuer  E  ^  ^  exactement. 


Nous  trouvons 


CO 


=  -  f.4  M  H  ? 


■R  (2) 


(?'M) 

4>M(0 


(69) 


Pour  evaluer  Eg'1-''  nous  pourrions  proceder  d'une  maniere 
semblable.  Cependant  la  sonme  est  suffisament  compliquee  pour  qu'un 
calcul  numerique  direct  sort  assez  fastidieux.  Gorame  nous  nous  attendons 
a  ce  que  E  soit  assez  petit  devant  nous  pouvons  tolerer  une 

grandeur  quolque  peu  approchee .  Nous  utilisons  done  un  rapprochement 
du  a  Ludwig.^ ^  Dans  cette  approximation,  lorsque  les  expressions  (51) 
pour  les  fonctions  de  Green  se  substituent  dans  (67),  les  frcjci  ),  q\il 
se  trouvent  dans  le  denominateur  sont  remplaces  par  leurs  grandeurs 
moyennes  prises  sur  le  spectre  des  frequences  du  cristal  pur.  On  peut 
alors  les  enlever  des  sornmes..  Les  sommes,  qui  fesultent  peuvent  etre 
evaluees  exactement.  Dans  ce  cas,  cette  approximation  est  equivalente 
au  remplacement  de  1' expression  exacte  pour  par  1' approximation 


G,fUr) 


(70) 


ou 


Ludwig,  J.  Pliys.  Chem.  Solids  9>  283  (1958)- 
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Ui  -  —  H  tfttii)  =  4  4="(r.) 

r  3N  U 


(71) 


Quoique  l'expression  donnee  par  (JO)  paraisse  grossiere  elle  donne  de  bons 


resultats.  Par  exemple,  si  l'on  s'en  sert  pour  e valuer  E 


(1) 


on 


trouve  le  mime  resultat  mais  le  coefficient  dans  (69)  est  1.5  au 
lieu  de  1.4448.  En  utilisant  (70)  et  (67),  nous  ob tenons 

e“’  =  ?  XZ  f  ,wv)  g  uun  f  «•) 

~  I 

2K 


«ji  v  8 

-  z:  t 

*  r  r 


(72) 


=  1  )  p"(0. 

Enfin,  de  (69)  et  (72)  nous  trouvons  que  l’energie  de  la  distortion 
statique  du  reseau  est 


Es  = 


.  |.  q  qq  9  iLiHl!  [l  -  O.Z53S  ]  .  (73) 

45"  (r.)  ■ 


II  nous  faut  maintenant  evaluer  la  contribution  dynamique  a 

✓ 

l'energie  libre  de  Helmholtz.; 
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La  somme  S^,  (65V.  devient  _ 

S,  =  12  N  4>"  «)  =  3NH  <oj 


Cependant-  que  la  somme  est  donnee  par 


Sz  =  2M  [  p"  (r. )  +  1  p'(r.)  ] 


Ainsi  la  contribution  dynamique  &  l'energie  litre  de  Helmholtz 
peut  s'ecrire  , 

=  3N  JcT  JL  0.92.51  f  Si  ^  St  1 
k.T  3  N  M  x 

.=  3  N  IcT  iw  -i-  0.92,51  [to}  4-  It  ]  1 

kT  ■  3NM 

s  +  4'  .  _ 


^  Z  MwJ 

Pour  notre  modele  la  grandeur  de 


<$>  4-  A  <&.  e» 


<£>  +  A  <§>t  =  &N  4=  (o +  ia.  pfi 


En  reunissant  les  differentes  contributions,  nous  pouvons 
ecrire  l'energie  libre  de  Helmholtz  du  cristal  deforme  dans  la  limite 
des  hautes  temperatues  sous  la  forme 

F  =  6N4><v.)  +  3NLTF  yjsi^  + 

IcT  (77) 

+  12  p(r.)  -  I.  HMH  ?  jfboF  r,  .  0.Z51S 

4>"M 

*  —■  “4r  2.M  Cp“f.)  +•  7-  p'(Oj. 
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Pour  1' analyse  4e  IV,  il  est  plus  convenable  d' exprimer  ce 
re suit at  sous  une  autre  forme.  Nous  supposons  que  le  deplacement 
instantane  Vji)  puisse  s'ecrire  comme  une  fraction  de  sa  grandeur 
a  I'equilibre  y^*'  (l  ) 


v«  U)  *  Xv.'"(i), 


(78) 


Ainai  l'energie  de  la  distortion  statique  peut  s’ecrire 


t  x1  r  v:'u)v’‘'(r)  *-■ 


JU< 

*  I1 


(79) 


+  X.  Z  IJlX'U)  +  fx  Z  v?U)V}"U’)  «■ 


JU1 

«/> 


D'apres  (53)  on  peut  recrire  cette  expression 

E  »  CX-  i  X1)  Z  , 


(80) 


et  nous  voyons  qu'elle  a  sa  grandeur  maximum  a  I  ,  comme  il 
convient . 

De  plus,  il  nous  faut  calculer  l'energie  de  la  distortion 
statique  dans  le  cristal  pur.  On  trouve 

Ej'1  =  I  x1  z  4>,f  (»•)  vJTU )  v?M .  (81) 


U' 


f 


On  pent  simplifier  lr evaluation  de  cette  sornme  si  nous  nous  servons  de 


Es°5t  =  -  i  y-{Z  +  } 

4  fl1  r  1 

=  A Z  I.HMHS  I  f  |  _  0,151$  ]  -(82) 

4>"  (rt )  L 

XL' rvcuvf'u)  v("u) 

ou  nous  a, vo ns  iitilise  (69).  Afin  d'evaluer  (82)  notis  avons  utilise 
1 '  approximation  contenue  dans  (70)  •  Comme  le  calcul  est  seniblalle 
a  celui  que  nous  avons  deja  fait  pour  evaluer  11  n'est  pas 

Jd 

necessaire  de  reproduire  les  details  ici.  Le  results, t  est 

Z> <fU)v?u)v?(U  *  JL  f«(r.) 

f  .  r  1  -  -  tlsm  j 1  (83) 

*■  H  <£"(r ,)  J 


En  se  servant  de  (82)  et  (83),  nous  trouvons  enfin 

Es  =  A1  l.HHH  ?  fi~  Q.LS1S  ?"(r.)  1 

....  III  >  4  ^ 


4>"M 


-  r  ± 


p"(re)  (  p^nzn-  1  i2- 

1  4>"(o'  H  4>>.)  J 
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ilombre  de 

r~™ - - 

<j\A  0  1 

3  .69  .26  • 

5  -54  .3^2 

7  . 42  . 378 

10  . 282  . 376 


Tableau  I 


i  comine  premiers  voisins 
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2 

3 

4 

5 

.044  ’ 

.005 

.001  ■ 

— 

.009 

.017 

.002 

— 

•  157 

•  039 

.006 

-- 

•  23 

.085 

.021 

.00  6 

ohm  cm 


ohm -cm 
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